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ABSTRACT 


Motivated  by  the  observation  that  the  photon  asymmetry  in 
radiative  muon  capture  is  determined  by  terms  of  0  ( 1/m 2)  in 
an  expansion  in  powers  of  the  nucleon  mass  m,  a  new 
calculation  of  radiative  muon  capture  in  40Ca  was  made  which 
is  consistent,  within  the  standard  theory,  to  0(1/m2)  -  The 
Hamiltonian  was  expanded  via  a  Foldy-W outhuy sen  reduction 
through  0(1/m2)  instead  of  the  usual  0(1/m)  to  allow  for  the 
first  time  a  calculation  of  the  0 ( 1 ) xO  (1/m2)  contributions 
to  the  rate.  The  results  show  that  the  Q(1/m2)  terms  are 
definitely  necessary  but  that  the  most  significant  ones 
arise  from  the  0  (1/m) xO  (1/m)  contributions  which  can  be 
obtained  from  the  usual  0  (1/m)  Hamiltonian. 

A  number  of  nuclear  effects  were  also  considered 
including  an  improved  giant  dipole  resonance  model  and  a 
model  using  Hartree-Fock  wave  functions  including  spin -or bit 
coupling.  As  in  most  previous  calculations  the  ratio  of 
radiative  to  ordinary  rates  and  the  photon  asymmetry  are 
relatively  insensitive  to  the  details  of  the  model,  though 
the  absolute  rates  can  be  quite  sensitive.  Fits  to  the 
available  data  to  extract  the  induced  pseudoscalar  coupling  gp 

were  also  done. 

In  an  improvement  to  the  closure  approximation,  the 
photon  spectrum  was  expanded  about  the  average  maximum 
photon  energy  km  and  the  correction  terms  evaluated  using 
for  one  a  modified  Thomas— Reiche—K uhn  sum  rule.  The 
resulting  rate  is  much  less  dependent:  on  km  than  the  usual 


v 


closure  result-  The  ratio  km/V  appropriate  for  closure 
calculations,  with  V  the  average  neutrino  energy,  was 
determined  and  found  to  be  approximately  constant  and, 
correction  terms  included,  somewhat  higher  than  values 
previously  used-  By  similar  technigues  a  consistency 
relation  was  derived  which  can  be  solved  to  explicitly 
estimate  nphysicaln  values  of  and  V  - 
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I. 

INTRODUCTION 

Eadiative  muon  capture  (EMC)  is  a  relatively  rare  process 
first  observed  in  iron  by  Conforto  et  al  (Co  62).  The 
measured  branching  ratio  for  radiative  capture  of  roughly 
10”4  times  the  ordinary  capture  rate  was  later  verified  by 
Chu  et  al  (Ch  65) ,  using  a  copper  target.  Subsequent 
experimental  effort  largely  focussed  on  the  40Ca  nucleus. 

The  40Ca  photon  spectrum  was  initially  measured  by  Conversi 
et  al  (Co  64)  at  CEEN ;  the  photon  spectrum  plus  asymmetry 
have  since  been  obtained  by  Di  Leila,  Rosenstein  et  al 
(Di  71,  Eo  73)  at  Nevis,  and  Hart  et  al  (Ha  77)  at  SEEL.  An 
experiment  proposed  for  TEIUME  (Ha  75)  will  look  at  the 
photon  spectrum  and  asymmetry  for  40Ca  and  20sPb. 

Early  theoretical  studies  (Hu  57,  Ma  59,  Op  64)  were 
confined  to  radiative  capture  on  a  free  proton  since  no 
complicating  nuclear  structure  effects  appear.  Experiments 
using  lh  targets  are  barely  feasible  because  of  the 
extremely  low  absolute  capture  rates,  so  cheoretical  efxorts 
naturally  shifted  to  complex  nuclei.  One  of  the  first 
theoretical  studies  of  radiative  capture  on  nuclei  was 
carried  out  by  Primakoff  (Pr  59)  ,  who  used  the  closure 
approximation  and  considered  only  the  muon  radiating 
diagram.  Eood  and  Tolhoek  (Eo  65)  made  a  detailed 
investigation  of  EMC  in  160  and  40Ca  using  the  impulse 
approximation.  Eeferences  to  the  "standard"  theory  of 
radiative  muon  capture  made  here  really  refer  to  the  Eood 
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and  Tolhoek  work.  Fearing  (Fe  66)  applied  Foldy  and 
Walecka's  (Fo  64)  giant  dipole  resonance  model  to  EMC,  and 
more  recently  Rood  et  al  (Eo  74)  considered  the  effects  of 
including  Coulomt  and  strong  nuclear  potentials.  The  advent 
cf  the  “meson  factories"  has  opened  up  the  possibility  of 
studying  EMC  on  very  light  nuclei,  for  which  newer 
calculations  also  exist  (Hw  78A,  Hw  78B). 

Before  this  investigation  was  started  and  prior  to  the 
experiment  of  Hart  et  al  (Ha  77) ,  a  serious  discrepancy  had 
existed  between  the  experimentally  measured  and 
theoretically  predicted  values  of  the  photon  asymmetry  for 
40ca .  Di  Leila  et  al  (Di  71)  had  determined  the  asymmetry  to 
he  0C<- 0.32±0.48,  averaged  over  the  interval  57<k<75  Mev. 

This  limit  on  OC  was  at  variance  with  the  standard  theory 
prediction  ^=+0. 75,  averaged  over  the  appropriate  photon 
energy  interval.  The  sign  discrepancy  has  now  been  largely 
resolved  by  the  recent  measurement  of  Hart  et  al  (Ha  77) 
yielding  0C=  +0. 90±0. 50,  though  there  is  still  room  for  more 
precise  numerical  values. 

In  1S7 5,  Fearing  (Fe  75)  proved  that  in  the  standard 
theory,  assuming  g5=0,  the  entire  first  order  correction  to 
OC  (in  zeroth  order,  ^  =  1)  comes  from  0  (1/m2)  terms  in  the 
square  of  the  transition  matrix  element.  Since  in  previous 
calculations  these  0(1/m2)  terms  had  not  been  retained 
consistently,  the  present  study  was  undertaken  with  the 
primary  objective  of  calculating  the  EMC  photon  spectrum  and 
asymmetry  consistently  to  0  (1/m2) .  At  the  same  time,  su^h 
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theoretical  improvements  and  supplementary  investigations  as 
proved  feasible  were  also  undertaken  in  an  effort  to  explain 
the  above  discrepancy-  Thus  the  standard  theory  calculation 
was  carefully  checked  and  extended  to  include  all  0  (1/m2) 
terms  in  the  square  of  the  matrix  element/  the  giant  dipole 
resonance  model  calculation  of  fearing  was  re-worked  and 
updated,  and  a  modification  of  the  usual  closure 
approximation  involving  nuclear  sum  rules  was  applied  to  HMC 
in  order  to  gain  some  insight  into  the  average  maximum 
photon  energy  parameter  kw- 

Ihfc  calculation  was  done  in  impulse  approximation, 
starting  with  the  usual  four-fermion  interaction  Hamiltonian 
describing  ordinary  muon  capture  on  a  free  proton.  The 
radiation  field  was  minimally  coupled  to  the  hadron 
covariant  portion  of  the  Hamiltonian,  which  was  reduced  to 
non-rela ti vistic  form  via  the  Foldy- Wouthuysen  procedure, 
keeping  all  0(1/m2)  terms.  Transition  amplitudes 
corresponding  to  a  set  of  Feynman  diagrams  were  ontained  by 
pairing  parts  of  the  non-relativistically  reduced  hadron 
covariant  with  appropriate  lepton  covariants.  The  lepton 
covariants  were  reduced  to  non-rela tivi stic  form  bj 
neglecting  the  small  components  of  the  muon  wave  function. 
The  resulting  transition  amplitudes  were  squared 
analytically  and  finally,  upon  choosing  nuclear  wave 
functions,  observables  were  computed. 

The  detailed  presentation  of  our  work  parallels  the 
steps  in  the  method  of  calculation  just  outlined.  Chapter  II 
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describes  the  relativistic  Hamiltonian  and  its  reduction  to 
non-reiativistic  form-  The  lepton  covariants  required  to 
yield  amplitudes  equivalent  to  those  obtained  from  standard 
theory  Feynman  diagrams  are  enumerated.  Chapter  III  deals 
with  the  problem  of  squaring  the  total  transition  amplitude. 
In  several  phases  of  this  operation,  symmetries  were  invoked 
to  keep  the  number  of  terms  at  a  manageable  level.  In 
Chapter  IV  the  nuclear  physics  aspect  of  SMC  is  addressed. 
F.eviews  of  the  closure-harmonic  oscillator  and  giant  dipole 
resonance  models  are  combined  with  descriptions  of  our 
improvements  to  them.  In  particular,  Hartree— Fock  wave 
functions  are  introduced  and  recent  photoabsorption  data 
utilized  in  the  giant  dipole  resonance  model.  Chapter  V 
contains  our  results  for  the  4°Ca  photon  spectrum  and 
asymmetry  for  a  variety  of  nuclear  models.  Chapter  VI  is 
devoted  to  describing  a  modified  closure  calculation  of  SMC. 
This  approach  involves  nuclear  sum  rules  and  offers  a 
possible  means  of  determining  the  average  maximum  photon 
energy  parameter  k^-  Finally  our  conclusions  are  presented 
in  Chapter  VII,  along  with  a  few  observations  and 
suggestions  regarding  further  work.  Solely  technical  matters 
and  discussions  cf  peripheral  issues  have  been  relegated  to 
the  appendices. 


II. 


EfFECTIVE  HA  MILTONIAN  FCR  RADIATIVE  MOON  CAPTURE  BY  A 

NOCLEOS 


A.  Relativistic  Hamiltonian 

The  impulse  approximation  was  used  to  describe  nuclear 
radiative  muon  capture  (EMC) -  This  approach  excludes 
explicit  nuclear  mesonic  degrees  of  freedom  from 
consideration,  the  nucleus  being  regarded  as  a  collection  of 
nucleons  which  participate  individually  in  the  RMC  process, 
hence  to  calculate  observables,  such  as  the  photon  spectrum 
and  asymmetry,  it  is  first  necessary  to  calculate  the 
transition  operator  for  radiative  muon  capture  by  a  single 
proton.  This  was  done  by  making  a  non- re la ti vis tic  reduction 
of  the  relativistic  Hamiltonian  for  radiative  muon  capture 
ty  a  proton  to  obtain  an  effective  Hamiltonian  for  this 
elementary  process.  The  effective  Hamiltonian  for  the 
nuclear  case  was  subseguently  obtained  by  summing  the 
elementary  effective  Hamiltonian  over  all  protons  in  the 
parent  nucleus.  Correction  terms  of  order  1/A  (A  the  atomic 
number)  arising  from  nuclear  center  of  mass  motion  (Kr  74) 
are  expected  to  be  small  for  40Ca  and  were  omitted. 

The  Hamiltonian  for  a  proton  subject  to  strong, 
electromagnetic,  and  weak  forces  was  taken  to  be  :  1 

i  Notation  and  conventions  are  those  of  J.  D.  Bjorken  and 
S.D.  Drell,  Relativistic  Quantum  Mechanics  (McGraw  Hill  Book 
Company,  Inc.,  New  York,  1964). 
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H  -  Ho  +  He  t  Hw 


2-1 


where  : 


Ho  -  Y»  in  +  Yo  ¥•  F*  +  Vo  +  Vo  V 

He  -  cJp  (Ao-Yo  7*  a)  +  i  (Xn7n  +Xp7p)  Yo  cfrv 

Hw  -  ^  7"-  ($v  Y  ~F  7°^^  /■  ^  Y  Ys 


77}*  r  2m 


2-2 


The  strong  interaction  appears  in  He  as  V©t  &V,  a 
combination  scalar  potential  and  one  transforming  as  the 
zeroth  component  of  a  4-vector,  m  is  the  nucleon  mass,  and  ? 
is  the  momentum  operator,  Ihe  form  of  the  weak  Hamiltonian 
is  the  usual  one  evolved  from  the  current  x  current 
theory  of  weak  interactions,  G  is  the  Fermi  constant,  m*  is 
the  muon  mass,  and  g„  ,  gA  ,  gP ,  g5  /  gr  are  the  weak 

couplings.  The  isospin  operator  7-  changes  a  proton  into  a 
neutron.  g^=ntr-p^is  the  4-momentum  transferred  to  the 
hadronic  current  in  the  weak  interaction.  The  leptonic  current 
oZoc  includes  all  interactions  not  directly  coupled  to 
the  hadronic  current,  and  its  explicit  form  depends  upon 
which  particle  radiates.  Detailed  expressions  for  <£oc  are 
presented  later  in  the  chapter  when  we  make  the  connection 
between  our  Hamiltonian  formalism  and  the  set  of  Feynman 
diagrams  calculated  in  the  standard  theory.  The 
elec nromagnetic  Hamiltonian  H^-  contains  Coulomb  and 


radiation  fields,  including  those  due  to  the  anomalous 
magnetic  moments  of  proton  and  neutron.  The  electronic 
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charge  is  e,  A0  is  the  Coulomb  field  of  the  nucleus,  and  A 
is  the  photon  field*  X^=1.79,  X'/V=-1.91  are  the  anomalous 
magnetic  moments  of  proton  and  neutron  in  units  of  nuclear 
magnetons* 

F^  are  the  electromagnetic  field  strengths-  The  isospin  operators 
Xn  select  protons  and  neutrons,  respectively. 

E-  Hon-fielati vistic  Reduction 

The  fully  relativistic  Hamiltonian  H  was  reduced  to  a 
non-relativistic  Hamiltonian  using  the  Foldy-Wouthuysen 
transformation  (Fo  50,  Bj  64).  The  same  non-relativistic 
Hamiltonian  can  be  obtained  from  a  Pauli  reduction  but,  as 
Friar  has  pointed  out  (Fr  75)  one  must  be  careful  to 
preserve  the  overall  normalization  of  the  transition  rate 
and  deal  with  initial  and  final  states  in  a  symmetric  way  to 
preserve  Hermit icity-  When  working  to  0(1/m2) ,  as  here,  or 
to  higher  orders,  the  Foldy-Wouthuysen  procedure  is 
preferable  because  algebraic  complications  are  at  a  minimum. 
Furthermore  the  proper  overall  normalization  is  assured 
since  to  0(1/m2)  the  difference  in  normalization  of 
relativistic  and  non-relativistic  wave  functions  (Bj  64)  is 
exactly  cancelled  by  the  difference  in  relativistic  and 
non-relativistic  phase  space  factors  (Fe  78).  The 
Foldy-Wouthuysen  transformation  decouples  the  relativistic 
Hamiltonian  into  two  two-component  Hamiltonians:  one 
describes  nucleon  states  and  the  other  describes  antinucleon 
states.  Since  the  transformation  reduces  the  Hamiltonian  to 
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only  so-called  "even”  operators  which  do  not  mix  upper  and 
lower  components,  an  initial  nucleon  state  is  coupled  only 
to  a  final  nucleon  state-  The  influence  of  the  antinucleon 
states  is  manifest  in  the  complicated  form  of  the  effective 
Hamiltonian. 

The  initial  step  in  the  reduction  involved  decomposing 
H  into  odd  and  even  operators,  which  respectively  do  and  do 
not  mix  upper  and  lower  components  of  the  nucleon  wave 
function-  Thus  : 

&  & 

H  s  771  -f-  £o  -h  Cw  +  fro  +  frw 
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where  : 

£.  =  eft  A.  +  Yc  Vi-  K  +  (x*  ft  iXrfr)  Y.  1 tVF*/ 

<%  =  Y.Y-(P-eftA) -i-Zifaft+XrTACY-CVA^+Y-A) 

=  (Ygz  <r*fp  +  ^  f*+J*  y*  M  X* 

frw  -  T-  (fv  Xs) 

Vz  J  77Lc  •  2m  ' 

^  2-4 


To  enable  a  consistent  calculation  of  EMC  observables  to 
appropriate  orders  in  the  expansion  parameter  1/m,  H  was 
transformed  so  that  the  odd  operators  appeared  in  lowest 
order  as  0  ( 1/m 3)  terms  : 

//"=  v.(m+£)  +  s-J-tfafaefl-jL. 

2771  8771  8771 


2-  c 


The  explicit  form  of  the  non-rela tivi sti c  Hamiltonian 
was  determined  by  a  combined  computer  and  hand  evaluation.  A 
FORTRAN  program  performed  the  initial  operator  algebra. 


which  amounted  to  commuting  operators  into  a  standard  form. 
Each  vector  space  (isospin,  Dirac,  etc.)  having  elements  in 
the  expression  for  H”  was  represented  by  a  linked  list 
structure  composed  of  those  elements.  The  operator 
commutation  rules  appropriate  to  each  vector  space  were 
coded  and  then  utilized  as  the  commutators  were  evaluated. 
The  basic  algorithm  required  the  product  of  two  terms  to  be 
put  into  a  predetermined  canonical  form  which  subsequently 
was  either  printed  out  or  used  in  further  evaluations.  The 
result  was  an  expression  for  H"  involving  Dirac  matrices 
which  were  expressed  in  terms  of  (two-component)  Pauli 
matrices.  The  final  simplifications  and  projection  of  the 
nucleon  states  (upper  components)  were  done  by  hand. 

Throughout  the  procedure  outlined  above,  the  exact 
expressions  for  the  lepton  covariants  were  left 
unspecified.  This  generality  was  in  keeping  with  the  as  yet 
quite  general  form  for  the  effective  Hamiltonian,  which 
describes  several  different  weak,  electromagnetic,  and 
strong  processes.  The  next  task  was  to  isolate  those  parts 
of  the  effective  Hamiltonian  which  specifically  describe 
radiative  muon  capture. 

C.  Effective  Haailtonian 

To  begin  constructing  an  effective  interaction 
Hamiltonian  for  radiative  muon  capture,  the  non- relati visti 
Hamiltonian  H"  was  written  as  : 
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where  H”  is  the  kinetic  energy  plus  strong  potential  part  of 
H"  and  IL"  is  composed  of  terms  involving  the  weak  and 
electromagnetic  interactions.  The  states.  Green’s  functions, 
etc.  corresponding  to  H”  are  labelled  with  the  subscript  o. 
The  calculation  cf  quantities  corresponding  to  the  full 
Hamiltonian  H"  then  proceeded  using  the  quantities 
corresponding  to  H^’  as  a  starting  point  (Ne  66)  .  Thus  the 
full  nucleon  Green's  function  is  : 

&*(£)  =  &!(£)  Hz'  &*(E)  +  &*(£) 

=  &Ue)  +  &*(£)  Hz  '&*(£) 

2-7 

The  superscripts  plus  and  minus  denote  retarded  and  advanced 
Green's  functions,  respectively,  lor  the  transition  from  an 
initial  state  i  to  a  final  state  f,  the  full  T-matrix  is  : 

Tf .  CE)~  Ta(E)  H  (  <H>'\e,-{),  [ Hz  +  Hz  &*<E)  Hz]  <H,M(E,i)) 

The  photon  emission  and  weak  interaction  processes  appearing 

in  H"  were  calculated  to  first  order  in  perturbation  theory. 

This  is  a  good  approximation  since  the  (dimensionless) 

expansion  parameters  e2/hc  and  : 

o]  _  M  .  -L  .  J-  /77luC\z  &  X  6.5  X  10  "  (mwC2)2' 

4ir  Key  h  ) 

with  m^c2  the  intermediate  vector  boson  mass  in  Mev  and  G 
the  Fermi  constant,  are  both  <  0.0  1.  Thus  G+(E)  was  replaced 
by  G*  (E)  and  only  those  terms  in  T^.  in  which  e  and  G  each 
appeared  just  once  were  retained.  The  leading  term  To ^ 


represents  scattering  by  the  nuclear  potential  and  can  be 
dropped.  The  radiative  muon  capture  T-matrix  thus  takes  the 
form  : 


n 
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It  is  completely  specified  when  the  lepton  covariants 
appearing  in  HJJ  are  specified. 

The  initial  effort  was  directed  towards  selecting  the 
Hamiltonian  corresponding  to  the  Feynman  diagrams  shown  in 
fig.  1.  The  calculation  of  these  diagrams  comprises  the 
standard  theory  of  EMC  {Eo  65),  so  guite  naturally  was 
chosen  as  the  starting  point  for  our  investigation  of  higher 
order  terms  in  the  effective  Hamiltonian.  In  this  approach 
the  nuclear  potential  is  omitted,  hence  the  propagator  G+  (E) 
appearing  in  Eg.  2-9  becomes  a  free  particle  propagator. 

Each  of  the  £K  described  in  the  following  list  was 
associated  with  a  corresponding  Feynman  diagram  as  follows  : 
1.  incoming  muon,  weak  vertex,  outgoing  neutrino.  This 

covariant  was  associated  with  a  nucleon  radiating,  fig. 

1  (b)  and  1(c)- 

X*  (x)  =*  %  (x)  y<t  tt-Ys)  %  Cx) 

=  exp  ' x)  ~Mv  y<x  (/~  y$)  2-10 


2.  incoming  muon  which  radiates,  weak  vertex,  outgoing 
neutrino.  This  covariant  was  associated  with  the 
non-radiating  part  of  the  nucleon  current,  fig.  1(a). 
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Fig.  1.  Diagrams  contributing  to  radiative  muon  capture  in 
the  standard  theory.  Diagrams  (b)  and  (c)  include 
radiation  from  the  anomalous  magnetic  moment  of  the 
proton  and  neutron. 
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it* (x)  =J i4y  %  (x)  Yx  i/-Ys)i SF (x-y)  it  #(y)  fy-iy) 

=  -  4-  exp(i(wk-/x)-x)  jit  Yx +7ru) & liu. 

M  r  y-k)z- V  2-11 


incoming  maor,  weak  vertex  coupling  to  the  nucleons  via 
a  pion  which  radiates,  outgoing  neutrino.  This  covariant 
was  associated  with  the  non-radiating  part  of  the 
nucleon  current,  fig-  1{d). 


Z* (x)  =  -Zc77u  ff txpfa(k-z)'X)+i (y+z-it)' y)  7/?  G*(/+fc)7U 

/2k  JJ  J(Ztt)4  J//  Zz-  771-/ 


*  -Z67TW  6X0  U(V+k-/x)’X)  1/9  Zee  C/+Ys)' 
J2T  (jUrV)z-  771k 
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Figure  1  (e)  comes  from  the  requirement  of  gauge  invariance, 

->oc  -*  cc  x 

ie.  from  the  substitution  p  -*•  p  -e  A  in  in  Eg.  2-2. 
Mplitudes  for  the  Feynman  diagrams  in  fig.  1  were 
constructed  by  substituting  the  above  covariants  into  Eg. 
2-9,  which  represents  the  effective  non-rela tivistic 
transition  matrix  element  obtained  from  the  Hamiltonian  in 
Eg-  2-4.  In  simplifying  the  lepton  covariants,  it  was 
assumed  that  the  muon  was  at  rest  and  that  the  muon  wave 
function  was  constant  over  the  nuclear  volume.  The 
prescription  described  by  Luyten  et  al  (Lu  63)  was  used  to 
obtain  an  appropriate  average  value  for  the  muon  wave 
function- 

specification  of  the  lepten  covariants  completed  the 
construction  of  an  0(1/m2)  effective  Hamiltonian  which 
describes  radiative  muen  capture  by  a  proton.  Performing  a 
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sum  over  the  constituent  protons  of  the  capturing  nucleus 
then  yielded  a  Hamiltonian  for  nuclear  radiative  muon 
capture  of  the  form  : 

e*=  %Uh?,-v)?7l 

with  (Tl  the  lepton  spin  matrices  and  the  nuclear  operators 
contained  in  : 


77[  =  TTJ,  +  Cl  ■  7T£z 


2-14 


lo  facilitate  manipulation  by  the  algebraic  squaring  routines. 


77£  was  arranged  into  the  form  : 

m=Z T-G){£ (EitD^C^Fi  +Z(En-Di+Ci)<n.-6i },  <5 
L  /=/  J  c  i=/  J.=N+!  J 


-/  S-A 


7 


2-15 


where  F.  and  3,*  are  functions  of  the  vectors  £,  k,  /,  P/m,  and 

jL 

p  .  The  coefficients  ,  and  C;  correspond  to  terms 

cf  order  1,  1/m,  and  1/m2  respectively*  The  complete 

Hamiltonian  appears  in  Appendix  B.  The  connection  with  the 
work  or  Food  and  Tolhoek  (Bo  65)  is  made  by  taking  the 
contact  approximation  for  the  nucleon  propagator  Ge  (E) ,  ie. 
letting  the  mass  of  the  propagating  nucleon  go  to  infinity, 
when  this  is  done,  agreement  is  found  with  their  effective 
Hamiltonian  which  contains  only  a  selected  few  0(1/m2) 
terms.  Other  parts  of  H^'  are  found  to  be  in  agreement  with 
previous  Foldy- Wouthuysen  reductions  of  weak  (Fr  66,  Oh  66) 
and  electromagnetic  (Me  62)  Hamiltonians.  2 


2  The  last  term  in  the  expression  for  in  Eg.  (4)  of 
(Fr  66)  is  dimensionally  incorrect.  The  factor  ( U)  -  k2/2m) 
should  read  (cu/k  -  k/2m)  - 
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An  extension  of  the  standard  theory  which  includes 
Coulomb  and  strong  nuclear  potentials  has  been  considered  by 
Eood  et  al  (Eo  74)  ,  who  use  the  Martin  and  Glauber  formalism 
(Ma  58) .  Much  of  their  work  involves  wholly  numerical 
evaluations  of  EMC  amplitudes  which  are  checked  against  the 
results  of  Rood  and  Tolhoek  (Ro  65)  in  the  limit  of 
vanishing  nuclear  potentials-  The  original  intent  of  our 
work  was  to  calculate,  in  an  analytic  fashion,  higher  order 
corrections  to  the  photon  spectrum  and  asymmetry,  and  this 
was  accomplished  for  the  standard  theory.  The  task  of 
modifying  our  existing  calculation  to  include  Coulomb  and 
strong  nuclear  potentials  appeared  to  require  considerable 
further  effort,  however,  and  was  not  undertaken-  Some 
further  discussion  of  the  role  of  exterior  potentials  in  EMC 
appears  in  the  presentation  of  results  in  Chapter  V. 


III. 


THE  TBANSITION  MATEIX  ELEMENT  SQUABED 

A.  Squaring  Strategy 

The  calculation  of  observables  for  radiative  muon 
capture  requires  that  the  transition  matrix  element  be 
squared  and  summed  on  appropriate  spins.  The  steps  in  the 
derivation  are  the  usual  ones  encountered  in  squaring 
feynman  amplitudes  (B j  64) .  With  the  inclusion  of  phase 
space  and  kinematic  factors,  the  probability  for  radiative 
muon  capture  corresponding  to  the  nuclear  transition  la>-^ 
jb>,  the  photon  and  neutrino  being  emitted  with  momenta  k 
and  V,  while  the  polarizations  of  the  initial  muon  and 
photon  are  specified  by  3  and  X  is  (Bo  65)  : 

Fib  ( k  v,  1,  s)  dkd&k  d/2?  -  FFL  / fylav  Z k  (km,- A)  /Mu!  dAd&k  da? 

4(2*}  2t> 

An  average  value  1^1^  of  the  muon  wave  function  over  the 
nuclear  volume  has  been  extracted  from  the  matrix  element 
=<b|H6^|a>.  £  =  1/13  7  is  the  fine  structure  constant, 
rractical  difficulties  encountered  in  the  present  calulation 
of  1M/  |2  arose  as  a  result  of  the  large  number  of  terms 
considered,  which  are  sufficiently  numerous  to  warrant  being 
manipulated  by  a  computer  code,  for  this  reason  the 
algebraic  language  REDUCE2  (He  73)  was  used  to  do  a  large 
part  of  the  algebra.  In  the  current  chapter  the  methodology 
employed  in  obtaining  the  squared  matrix  element  is 
described  in  detail. 
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E.  Lepton  Spin  Algebra 

Eeginning  with  the  effective  Hamiltonian  : 


H  =  a-0l-Y)77[/.r 


3-2 


the  square  of  the  matrix  element  was  constructed,  taking 
account  of  the  initial  muon  polarization  by  introducing  the 
spin  projection  operator  (1  +  5*££)/2  to  the  immediate  left  of 
the  muon  spinor  : 

Z  /Mljz=  J.  Z  //.?(!-  <Z-  v)  m  (/+  CL-  S)XA  J* 

4  spins  /  (y 


spins 
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Ihe  most  general  scalar  form  for  77£  in  terms  of  the  lepton 
spin  matrices  (Tl  is  : 


;y  =  77£,  +  77[z 


3-4 


Substituting  this  form  in  Eg.  3-3  and  doing  the  trace 


algebra  yields  : 

Z  /MiJz=zlO-v-s)wrrb+(l+v-s)77irnizf-i(?+s)-7%zx7y2.f 

-Z/Re[77[, (r-i+*vxi)-77iz  -t-v-fy  S-^/J]  3-5 


spins 


Ihe  expressions  77f,  and  are  given  explicitly  in  Appendix  B. 
Eg.  3-5  was  evaluated  analytically  with  the  help  of  the 
computer  language  REDUCE2 .  Before  describing  the  explicit 
evaluation  though,  some  general  si  mplif  ica  tions  of  the 
squared  matrix  element  are  discussed. 
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C.  Symmetries  in  tSoe  Squared  Matrix  Element 

The  squared  matrix  element  exhibits  particular 
symmetries  which  can  be  utilized  to  simplify  it 
considerably.  The  manifestation  of  these  symmetries  may 
require  some  idealizing  assumptions,  and  where  this  is  the 
case,  specific  mention  is  made  of  these. 

Nucleon  Spin 

Consider  that  part  of  the  squared  matrix  element  which, 
in  closure  approximation,  has  the  form  : 

<nz  V(j)Z  &(*)%(*)/*> 

+  I 

where  0 ( j )  and  0  (k)  are  single  particle  operators 
independent  of  spin.  It  is  assumed  that 

1.  ja>  is  a  Slater  determinant  of  single  particle  wave 
functions  including  spin-orbit  coupling  and  that 

2.  1 a>  has  closed  subshells. 

A  straightforward  calculation  using  the  techniques  of 
angular  momentum  algebra  (Ed  57)  shows  that  the  one-body 
(j=k)  part  of  Ifyl 2  is  identically  zero  and  that  its  two-body 
(j*k)  part  is  proportional  to  the  product  of  3-j  symbols 

f  j,  a  i  ]  / A  A  J'W  J  A  i  \ 

which  also  vanishes.  Here  1;  and  lz  are  orbital  angular 
momenta  of  the  single— particle  states.  Hence  terms  in  our 
squared  matrix  element  involving  only  one  nucleon  O'  are 


identically  zero. 


In  the  case  where  closure  is  not  used,  Luyten  et  al 
(Lu  63)  find  that  in  addition  to  the  above  assumptions  one 


requires  no  spin-orbit  coupling  in  order  that  <1X0=0. 
Integration  on  Neutrino  Angle 

When  the  initial  muon  has  polarization  3  ,  the 
available  vectors  which  can  be  used  to  construct  observables 
after  integrating  over  neutrino  angle  and  summing  on  spins 
are  3  and  k.  As  a  result  the  differential  photon  spectrum 
exhibits  the  following  dependence  on  3  and  k  : 


3-7 


with  X  the  photon  asymmetry.  Hence  is  invariant  under  the 


transf ormation  «5,k-*-5,-k.  If  the  final  and  initial  nuclear 
states  have  definite  parity,  the  transformation  r-*-r 
applied  to  the  matrix  element  will  not  affect  the  wave 
functions  appearing  in  |  2  since  each  state  appears  twice, 

however  P  is  is  changed  to  -P.  Making  these  transformations 
on  |  J  2  and  now  taking  V-*-V  has  the  effect  of  making  each 
term  in  the  matrix  element  an  eigenstate  of  the 
transformation  operator  which  changes  V-+-V  and  has 
eigenvalue  +1  or  -1.  Thus  from  an  operational  point  of  view, 
all  terms  exhibiting  odd  behaviour  under  the  transformation 
?,k,3,P--»,-k,-5,-P  vanish  when  integrated  over  neutrino 


solid  angle.  Note  that  as  a  consequence  of  this,  the  final 
result  for  the  squared  matrix  element  integrated  over 
neutrino  angle  will  not  contain  any  vector  triple  products. 
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Ee lations  Between  Squared  Nuclear  Matr ix  Elements 

Define  the  following  set  of  squared  nuclear  matrix 
elements  for  ordinary  muon  capture  : 

/MalZ=  Z  (%f)  Nfy  /Ci/z  Oa(l)  exp  (-i  ?tl-Ai)/a.)>l 

r  a  4ir  i*,  r  /  3.8 

where  p  can  he  V,  A,  or  P  with  : 

Ov  U)  =  1M  ,  Oa  U)  =  -L  7U)  .  Or(i)°  V-  CTU) 

J  \f3  3-9 

Then  the  following  relations  between  matrix  elements 

(Eo  64)  : 

/MvlZ=  /MaI*‘  /Mr/* 

3-10 

hold  exactly  in  a  single  particle  model  without  spin-orbit 
coupling,  and  hold  to  a  good  approximation  in  cases  where 
the  spin  dependence  of  the  inter-nucleon  forces  is  fairly 
weak  (lu  63,  Fo  64,  Go  74).  Applied  to  radiative  muon 
capture  (Eo  65,  Fe  66) ,  the  relations  appear  as  : 

Z  (kt-kf ’/</>/£  7-  a)  exp  A)  Mi)  I  a)  /  * 

=  3  Z  (km-tif/ <b/Z  7-  a)  exp  (-J.  -A)  I  a}  1 2 

ib  1  7 

XL  (kn-kf A-  <&/Z  Z  U)  exp  (~i  Zt  A)  Mo /Z>  B*<h  /Z7-fi)exp(-iS<u.-A)crU)  A> 

=  Z  A- B’&Zt-k)2' KA/Z Z- d) exp(-£S,t  Ai)Ia.yr  3-11 

-  -»  -»  — ►  -*■ 

where  s ^  =  (k+p)^  •  and  A'  B  are  arbitrary  vectors.  Note  tnar 
the  net  effect  of  the  relations  3-10  is  to  remove  the 
nucleon  spin  operators  from  the  matrix  element- 


21 


Dm  Evaluation  of  the  Square  of  the  Matrix  Element 

With  the  help  of  the  simplifying  relations  detailed  in 
the  previous  section,  the  analytic  evaluation  of  Eg.  3-5  was 
carried  out.  The  algebraic  manipulations  were  done  in  three 
distinct  parts,  each  part  being  accomplished  by  a  BEDUCE2 
program  written  specifically  for  the  task. 

The  initial  phase  of  the  work  involved  classifying  the 
terms  in  the  effective  Hamiltonian  according  to  : 

1.  their  behaviour  under  the  transformation  V,k, 

-y,-k,-5,-P,  and 

2.  the  number  of  nucleon  (f's  appearing  in  each  term. 

All  possible  forms  emerging  from  this  classif ication  scheme 
(in  which  the  Coulomb  gauge  was  used  and  the  photon 


was  made  to  appear  explicitly) 


were  enumerated  and  substituted  into  Eg.  3-4.  A  list  of 
these  forms  can  be  found  in  Appendix  B.  Eg.  3-5  was  then 
evaluated  and  the  simplifications  cf  the  preceding  section 
associated  with  points  1.  and  2.  above  were  done.  Each 
product  of  forms  in  this  intermediate  expression  represented 
a  large  number  of  terms  in  the  sguared  matrix  element. 
Elimination  of  the  photon  polarization  vector  was  effected 
by  the  relation  : 


~  3-12 

-*•  (%) 

which  follows  directly  from  the  explicit  expression  for  € 


This  first  stage  of  the  sguaring  procedure  was  carried  out 
by  the  program  SQUARE1. 


22 


The  next  step  was  to  substitute  the  original  terms  back 
into  the  expression  for  the  squared  matrix  element ,  which  at 
this  stage  was  simplified  insofar  as  the  symmetry  properties 
of  the  forms  had  been  fully  exploited.  The  re-substitutions 
could  not  be  done  all  at  once  as  the  number  of  terms 
generated  by  the  substitutions  far  exceeded  REDUCE2's 
storage  capacity.  Rather  the  squared  matrix  element  was 
sub-divided  and  the  segments  processed  one  at  a  time.  The 
program  SQUARE2  handled  the  re-sub  stit utions  and  the 
elimination  of  nucleon  (T's  via  the  relations  in  Eqs.  3- 11a 
and  3-r  1  1b. 

The  final  analytic  form  of  the  squared  matrix  element 
was  produced  by  the  program  SQUARE3.  An  integration  over  the 
(experimentally  unobserved)  neutrino  solid  angle  was 

A  A 

performed,  introducing  y=k*v?  as  an  integration  variable.  The 
integrals  of  various  powers  of  y  multiplied  by  a  squared 
nuclear  matrix  element  were  systematically  labelled.  A 
feature  of  the  REDUCE2  language  permitted  output  of  the 
squared  matrix  element  from  SQUARE3  in  FORTRAN,  thus  greatly 
facilitating  further  numerical  work. 

In  our  final  expression  for  the  squared  matrix  element 
all  terms  of  0(1/m2)  and  the  0(1/m3)  terms  coming  from  the 
product  of  0  (1/m)  and  0(1/m2)  terms  in  the  effective 
Hamiltonian  were  retained.  These  included  the  so-called 
•’velocity  terms",  ie.  those  terms  involving  nucleon  momenta. 
Appendix  C  contains  an  explicit  expression  for  the  squared 
matrix  element  which  is  fully  specified  up  to  the  nuclear 
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matrix  elements.  A  number  of  nuclear  models  permitting 
evaluation  of  the  latter  quantities  is  presented  in  the  next 
chapter. 


IV 


NUCLEAB  MODELS  FOE  RADIATIVE  MUON  CAPTUBE 

A-  EMC  Models  for  Complex  Nuclei 

In  this  chapter  two  nuclear  models  for  radiative  muon 
capture  and  our  improvements  to  them  are  discussed.  The 
standard  closure-harmonic  oscillator  (CHO)  model  employed  by 
Food  and  Tolhoek  (Eo  65)  and  the  giant  dipole  resonance 
(GDE)  model  of  Foldy  and  kalecka  (fo  64)  as  applied  to 
radiative  muon  capture  by  Fearing  (Fe  66)  are  reviewed. 
Included  in  the  discussion  of  the  har monic-oscilla tor  model 
is  an  improvement  to  it  in  the  form  of  a  more  sophisticated 
shell  model  (SBL)  wave  function  due  to  Shao,  Eassichis  and 
Lomon  (Sh  72) .  In  addition,  a  variant  of  the 
closure-harmonic  oscillator  model  involving  nuclear  sum 
rules  is  mentioned  briefly.  A  cocre  complete  discussion  of 
this  latter  approach  as  applied  to  radiative  muon  capture  is 
contained  in  Chapter  VI-  The  standard  closure-harmonic 
oscillator  model  is  discussed  first. 

E-  Closure  -  Harnonic  Oscillator  Model 

Si m pie  Wave  Functions 

Recalling  Eg-  3-1,  which  gives  the  probability  for 
radiative  muon  capture  corresponding  to  the  nuclear 
transition  Ja>-*|b>  : 

Pailky^.l)  dA  dsi/k  d-ft-v  ^  /**£  A  (A*.- A)  /Mu!  dA  d-&A  ^  ^ 
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where  M^A=<h  |  He^|  a>#  it  is  apparent  that  the  exact 
calculation  of  : 

JL  (k™~k)Z /Mu I* 

Zb 

requires  a  knowledge  of  the  wave  functions  of  the  accessible 
states  in  the  daughter  nucleus*  for  40Ca,  the  maximum  energy 
available  to  excite  the  daughter  nucleus  4 °K  is  : 

Em  =77^'  Ebe  -  (m*-  772 p)  +  Ec  -  //0. 4 

4-1 

in  which  E^-  is  the  muon  binding  energy  and  Ec  is  the 
Coulomb  energy  difference  between  initial  and  final  nuclear 
states.  Since  it  is  expected  that  most  of  the  transition 
strength  lies  in  transitions  to  relatively  low-lying  states 
in  the  daughter  nucleus,  the  use  of  the  closure 
approximation  is  appropriate  (Pr  59).  Closure  involves 
replacing  the  sum  cn  all  possible  final  states  1 b>  in  Eg. 

3-1  by  a  sum  on  a  complete  set  cf  states  : 

Z  IhXbl  =  1 

&  4-2 

The  introduction  of  superfluous  states  E >E^  in  this  manner 
is  not  a  cause  for  concern  if  the  majority  of  transitions 
are  to  states  of  sufficiently  low  energy,  while  the 
advantage  gained  is  that  specific  final  state  dependent 
quantities  are  replaced  by  suitable  averages.  Thus  km  is 
replaced  by  an  average  maximum  phctcn  energy  km 
corresponding  to  an  average  exitation  energy  of  the  daughter 


nucleus  : 


Mm  —  77^u.  ~~  ( 771n  ~  771  p)  £&E  ( £.<l)  *  Mm 
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with  the  associated  notational  change  : 


( V+k)al>  V  +  k  ®  5 
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la ,  are  the  energies  of  the  initial  and  final  nuclear 
states,  respectively.  It  now  remains  to  evaluate  products  of 
nuclear  matrix  elements  of  the  general  form  : 

<o,y  <(h  >*  =  <a/Z  Jt  (p  T+  (k)  txp  (-/  s-  (Aj  -Ak ))  0,  (j)  ol(M)  J  a) 

/>  *  4-  5 


which,  since  they  are  ground-state  expectation  values  of 
two-body  operators,  can  be  calculated  relatively  easily  in  a 
nuclear  model.  The  sums  on  j  and  k  run  over  all  occupied 
nucleon  states  of  the  parent  nucleus.  The  operators  0/  and  0 
are  in  general  functions  of  nucleon  momenta  and  spins. 

Here  <  IX  1>  will  be  referred  to  as  the  main  nuclear  matrix 
element  squared  since  it  dominates  numerically.  For  40Ca, 
the  result  for  the  main  nuclear  matrix  element  squared  in 
the  closure-harmonic  oscillator  model  is  : 

z 

</></>*=  zo-  (ZO  t  2.5 if-  O.zsf  +  0.  03/25  f. )  &'  ^ 

4-6 


with  ^  2=b2  (s- s)  =  (sb)  2  and  b  the  oscillator  parameter. 
Explicit  details  of  the  evaluation  of  the  products  of  the 
nuclear  matrix  elements  required  here  can  be  found  in 
Appendix  D. 

Hartree-Fock  Wave  Functions 

A  straightforward  improvement  to  the  closure-harmonic 


oscillator  model  was  achieved  by  replacing  the  simple 
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single- particle  oscillator  wave  functions  for  the  nuclear 
matrix  elements  <1><1>  and  <1><P>  with  the  more  realistic 
ones  of  Shao  et  al  (Sh  72) ,  calculated  from  the 
Feshbach-Iomon  potential  using  a  Hartree-Fock  technique.  3 
Since  the  SB1  wave  function  includes  spin-orbit  coupling, 
the  relations  | My | 2  =  | M^ | 2  =  | 2  adapted  for  EMC  (Egs.  3-11) 
can  be  expected  to  hold  only  approximately.  An  explicit 
calculation  shows  |M>J12  =  |M/|2  to  better  than  3%  and 
I | 2= | j  2+4,  where  the  correction  A  involves  the  tensor 
coupling  between  M p  and  M*.  A  is  identically  zero  when  JM^j2 
is  integrated  on  photon  and  neutrino  solid  angle,  and  is 
estimated  to  be  very  small  otherwise.  These  results  tend  to 
indicate  that  the  approximation  |M^ |2= |M^ | 2= | Mp| 2  is  a  good 
one,  and  we  have  used  it  here. 

Shao  and  collaborators  have  expressed  their  wave 
functions  in  a  basis  of  harmonic  oscillator  functions  each 
having  the  same  oscillator  parameter,  thus  permitting  the 
main  nuclear  matrix  element  squared  to  be  calculated 
analytically.  Working  with  a  basis  of  oscillator  functions 

Hf 

meant  that  the  actual  form  of  the  result  for  <1X1>  was 

known,  so  that  the  problem  became  one  of  determining 

coeff icients.  This  task  was  straightforward  as  the  angular 

momentum  algebra  was  not  particularly  complicated,  and 

closed  form  expressions  exist  for  the  radial  integrals 

3  <1><p/  was  expressed  in  terms  of  the  two-body  part  of 
<1X1>*and  was  calculated  concurrently  with  it.  The  other 
nuclear  matrix  elements,  being  of  lesser  numerical 
significance,  were  calculated  using  the  simple  wave 
functions  described  in  the  previous  section.  See  Appendix  D 
for  details. 
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(De  66)  -  The  large  number  of  terms  to  be  multiplied  and 

summed  dictated  a  computer  evaluation.  Details  of  the 

calculation  are  to  be  found  in  Appendix  E.  The  result  : 

z 

</></>*=  10-  (zo -3. OSy-t  S./S3  if-  J.  (,Z7f+  0. 4/4f)e^ 
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may  not  be  compared  directly  with  the  result  for  simple  wave 
functions  given  in  Eg.  4-6  since  in  this  latter  equation 
b=2. 03  fm  is  appropriate,  while  in  Eg.  4-7  b=1,90  fm  should 
be  used.  Different  methods  for  fixing  b  have  resulted  in 
these  two  different  values.  Rood  et  al  (Ro  65)  determined 
b=2.03  fm.  by  regurring  that  the  simple  harmonic  oscillator 
wave  functions  reproduce  the  40Ca  r.m.s.  charge  radius 
obtained  from  electron  scattering  experiments.  Shao  et  al 
(Sh  72)  determined  b=1.90  fm.  by  requiring  that  their 
hartree-fock  wave  functions  yield  optimum  values  for  the  1 60 
r.m.s.  charge  radius  and  binding  energy- 

It  may  be  noted  that  the  expression  in  Eg.  4-7  contains 
a  ^2  term  whereas  the  expression  in  Eg.  4-6  does  not.  The 
explanation  for  this  difference  of  form  (Be  71)  is  to  be 
found  in  the  property  of  factorization  of  simple  harmonic 
oscillator  wave  functions  into  products  independently 
describing  relative  and  center  of  mass  motion,  combined  with 
the  fact  that  for  40Ca,  N=Z.  Explicitly,  if  rsospin  is  a 
good  quantum  number,  the  closure  result  for  the  main  nuclear 
matrix  element  sguared  can  be  expressed  : 


<7>  0>  * =  <«-l Z  Ot  (J)  exp  U  S-Aj)  Z  XU)  exp  (-J.  s£)/a.> 


29 


<7><7>*=  exp(i$fij)  £  Ji U)  up C-M  S-/Z)M)  4-8 

lor  nuclei  with  N=Z  the  tern  of  second  order  in  ^7=sb  in  Eg. 
4-8  is  : 

~  z  K&lX*  (l)  S'Ai  X'  (j )  S'Aj  / &y 

9  “  ~z  £  k&i S'  c r?~  r*)  / &y 

4-9 

where  E^  and  E^  are  center  of  mass  coordinates  of  the 
protons  and  neutrons.  For  the  simple  harmonic  oscillator, 
the  center  of  mass  motion  has  a  ground  state  wave  function 
proportional  to  : 

& xp  (~fj*  ( Rp-R" )  ) 

so  that  the  closure-harmonic  oscillator  model  gives  for  Eg. 

4-9  : 

_  J.  -t  * 

*■  z  ? 

*■ 

Since  the  general  form  of  the  result  for  <1X1  >  is  : 

</></>*=  Z~zf  t  +  C,  yZ+  Cz  w4+...  }  exn  (-  *l/z) 

1  1  J  4-10 

the  term  Z(-J/j?2)  is  entirely  accounted  for  by  the  expansion 

of  the  exponential  in  Eg.  4-10,  therefore  implying  C,  =0  for 

the  closure-harmonic  oscillator  model.  The  SBL  wave  function 

cn  the  other  hand  cannot  be  factored  into  terms  which 

individually  describe  relative  and  center  of  mass  motion,  so 

that  the  evaluation  of  Eg.  4-9  will  not  in  general  lead  to 

the  result  Z  (-f  ^  2)  #  and  hence  C,  will  be  non-zero. 

A  crucial  shortcoming  of  the  closure-harmonic 

oscillator  model  is  that  it  predicts  ordinary  capture  rates 

which  are  approximately  a  factor  of  two  larger  that  those 
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measured  for  160  and  40Ca.  A  significant  part  of  the 
discrepancy  may  be  understood  once  it  is  realized  that  most 
of  the  capture  occurs  through  the  first  forbidden  dipole 
matrix  elements  (li  49)  to  states  lying  in  the  region  of 
excitation  of  the  giant  dipole  resonance.  The  energies  of 
these  states  should  be  accurately  reproduced  by  the  nuclear 
model  because  the  capture  rate  depends  strongly  on  the 
maximum  photon  energy  (the  radiative  rate  goes  as  )  .  In 
the  simple  shell  model,  the  giant  dipole  resonance  does 
emerge  as  a  collective  feature,  but  at  too  low  an  energy,  so 
that  for  a  simple  shell  model  calculation  which  does  not 
employ  closure,  the  dipole  contribution  to  the  rate  will  be 
overestimated.  Presumably  the  closure-harmonic  oscillator 
model  rate  is  overestimated  for  the  same  reason,  A  means  of 
avoiding  this  difficulty  was  developed  by  Foldy  and  Waiecka 
(Fo  64)  in  an  investigation  of  ordinary  muon  capture.  Their 
approach,  subsequently  referred  to  as  the  giant  dipole 
resonance  model,  was  later  extended  to  radiative  muon 
capture  by  Fearing  (Fe  66),  and  it  is  in  this  latter  form 
that  the  model  is  reviewed  here. 

C.  Giant  Dipole  Eesonance  Model 

The  principal  idea  of  the  giant  dipole  resonance  model 
is  to  relate  the  dipole  part  of  the  muon  capture  rate  to  an 
integral  over  the  experimental  photo— absorption  cross 
section.  To  obtain  the  relation  it  is  first  assumed  that 
isospin  is  a  good  quantum  number.  Thus  for  a  nuclear  ground 
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state  j a>  having  total  isospin  T=0  (Fe  66)  : 


z 

tU? 

=  ~Z  Z  fe (fern~ fe)  /<Z/Z  TzU)  IXP  (~jL  5&l'Ac)/&yl  ^ 

&>'  ^3—  *  r 

~Z7T  Z  fetem’fef  Z  KZ/ZJzWu,  (SdAi)  vr^/*>r 

Xi*  77J*3  *  0*  4-11 

Ihe  sum  on  b*  runs  over  excited  T=1  states  of  the  parent 

nucleus.  Ihe  maximum  photon  energy  in  this  expression  is 
ab 

still  k^,  corresponding  to  the  proper  transition  energy. 


ab 


kjkmzti.  /<b!Z  j:  a)  txp  (-i  S*i  -Ai  )/<*>/ 


-rr>  3 


Hence  : 


ab 


fern  ~~  Zm  Zb  &. 
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with  Em  as  defined  in  Eg-  4-1. 

* 

In  their  examination  of  the  dipole  (1=1)  part  of  the 

ordinary  muon  capture  analogue  of  Eg-  4-11,  Foldy  and 

Walecka  (Fo  64)  noted  that  the  dipole  matrix  element  could 

be  written  as  the  product  of  the  matrix  element  of  an 

unretarded  part  S*bAi  and  an  elastic  form  factor  F^.  With 

3 

this  prescription  the  dipole  part  of  Eg.  4-11  becomes  : 

o&m-jrZ  k(k£-kf  lF.,lz(Sd.fkt!IZ  7xcaZ  M>/* 

6  W  77^3  *  4-13 

Ihe  closure-harmonic  oscillator  model  was  used  to  compute 
the  dipole  term  appearing  in  Eg.  4-11  and  the  approximation 
to  it  used  in  Eg.  4-13.  Besides  giving  a  check  on  the 
accuracy  of  the  approximation,  this  calculation  also  yielded 
a  check  on  the  formulae  used  by  Foldy  and  Walecka  (Fo  64) . 
Ihe  calculation  shows  : 
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z-rr  Z  z  K&'IZ  JzU)j,  {Sti/Ii)  y,  7j2a;)  IE)  1 2 

4Lo  771  4s  ff 

= i07i*{/-iyz+  M> -fa ->[**■  Two  -%*+...] txp(-n%)  4_14 


i  Z  /Fti!z  (Ssif  Kb'IZ  3k(i)Ju/A>r 
w 

JL 

=  lOrffl-iTf-t-so-rf-lio  Tf1,  b  3^00  ?1B }  &xp  (-l/z) 

Here  : 

y*-  ( sb)z s  bz[kL-ik(kn-k)(/-y)] 
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and  b  is  the  oscillator  parameter.  For  40Ca#  taking  b=2. 03 
fm  and  k  =90  Mev,  it  is  evident  that  Eg.  4-14  differs  from 
Eg.  4-15  by  at  most  a  few  percent  for  any  value  of  the 
photon  momentum  k,  so  that  Eg.  4-13  is  indeed  a  good 
approximation  to  the  dipole  part  of  Eg.  4-11. 

The  nuclear  matrix  elecent  which  appears  in  Eg.  4-13  is 
precisely  that  which  describes  nuclear  photo-absorption  in 
the  unretarded  dipole  approximation  : 

0}(£)=  i  7TZXZ  (£i-U)  Kb'/Z&U) %/*>/* £(£-Eb>-Ea) 

£6'  M  4-17 


"  -l—  /fa  (£/-c$)/Z fd£  f  &(£>  *  (£mr£-k)*  $(Em-£-td 

Z7T*X  /  <  £  TTJu.3 

X  [ (Em -  E) l-Zk  (£m-£-k)  (h y)j] 
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A  *  /Fu  (£tu>/  2  { a,  (k)  +  &  (k)  (/-  y) } 

4-18 

where  #(x)  is  the  unit  step  function,  ErtS  is  the  energy  of 
the  giant  dipole  resonance  and  the  f unctions  Q,  (k)  and  Qz  (k) 
are  defined  : 


Hence  the  dipole  part  of  the  capture  rate  is  related 
directly  to  the  total  photo-abscrption  cross  section, 
insofar  as  the  squared  nuclear  matrix  element  of  the  form  in 
Eg.  4-11  is  concerned. 

The  latter  qualification  is  necessary  in  that  generally 
the  nuclear  matrix  elements  are  functions  of  nucleon  spins 
and  momenta.  The  nucleon  spins,  however,  are  eliminated  from 
consideration  via  the  relations  between  squared  nuclear 
matrix  elements  presented  in  Sec.  Ill* C. 1,  so  that  the 
matrix  elements  appear  either  as  in  Eg*  4-11,  or  contain 
additional  factors  of  nuclecn  momenta  -  the  so-called 
"velocity  terms".  As  the  numerical  results  in  the  next 
chapter  will  show,  the  velocity  terms  play  a  lesser  role 
quantitatively  in  radiative  muon  capture.  In  the  giant 
dipole  resonance  model,  the  velocity  terras  and  the  other 
multipole  terms  (non-dipole  contributions)  were  evaluated 
via  the  closure- harmonic  oscillator  model,  using  either 
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simple  oscillator  or  SBL  wave  functions. 

D-  Modifications  to  the  Usual  Closure  Approximation 

In  discussing,  at  the  end  of  section  IV.  B,  the  failure 
of  the  closure- harmonic  oscillator  model  to  predict  the 
correct  dipole  transition  strength,  the  strong  km  dependence 
of  the  capture  rate  was  mentioned.  This  dependence  is  of 
immediate  concern  because  the  average  maximum  photon  energy 
A^is  not  determined  a  priori  by  the  model,  although  fixing  ^ m 
is  highly  desireable  from  the  point  of  view  of 
eliminating  uncertainties  due  to  the  nuclear  physics.  A 
similar  difficulty  arises  in  a  closure-harmonic  oscillator 
model  description  of  ordinary  muon  capture,  where  the 
parameter  to  be  fixed  is  the  average  neutrino  energy  .  A 
technique  employing  nuclear  sum  rules 

(Co  72,  Be  73,  Go  74,  Ko  76)  which  involves  modification  of 
the  simple  closure  approximation  has  been  successful  in 
largely  eliminating  the  dependence  of  the  ordinary  capture 
rate.  Recently  the  method  has  been  extended  to  radiative 
capture  by  Sloboda  and  Bearing  (Si  78)  ,  whose  work  appears 
as  Chapter  VI  here.  A  study  of  the  conditions  of 
applicability  of  the  modified  closure  approach  was 
undertaken  by  Rosenfelder  (Ro  77) ,  who  demonstrated  its 
suitability  for  the  muon  capture  problem. 

Ideally  a  combination  of  the  giant  dipole  resonance  and 
modified  closure- harmonic  oscillator  models  should  result  in 
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a  model  which  is  superior  to  either  of  the  two.  The 
realization  of  this  combination  meets  with  some  difficulties 
in  practice ,  however,  as  the  simple  sum  rule  arising  in  the 
modified  closure-harmonic  oscillator  model  does  not  emerge 
intact  in  attempts  to  modify  the  giant  dipole  resonance 
model.  The  proposed  combined  model  would  require  a  sum  rule 
for  the  other  multipoles  which  must  be  of  higher  order  than 
the  first  order  one  described  in  Chapter  VI.  This  is  so 
because  for  the  first  order  sum  rule,  the  dipole  part  of  the 
main  nuclear  matrix  element  squared  for  EMC  yields  the  same 
sum  as  the  full  main  nuclear  matrix  element  squared  itself. 
Higher  order  sum  rules  may  be  derived  (Er  70)  but  their 
complexity  makes  them  difficult  to  use. 

The  nuclear  models  outlined  above  complete  our 
theoretical  description  of  radiative  muon  capture.  In  the 
following  chapter  numerical  results  are  presented. 


V. 


BESDLTS 


A.  BMC  Observables 

By  way  of  introduction  to  the  numerical  results,  a 
summary  of  observables  for  radiative  muon  capture  follows. 

The  probability  for  radiative  muon  capture  from  an 
initial  nuclear  state  |a>  to  a  final  nuclear  state  |b>  was 
civen  in  Eg.  3-1  as  : 

P<d>  (ky  v?  3)  PAP/2&  Pft?  ~ xjPfar)4  ^ ^ ^ /Mu!  PA P-&&  P-fl-v 

^  3-1 

Specializing  to  doubly-closed  shell  nuclei  via  the  relations 
between  squared  nuclear  matrix  elements,  Eqs.  3-1  la  and 

3- 1 1b,  and  using  the  closure  approximation  detailed  in  Egs. 

4- 2  to  4-4,  the  total  transitior  probability  from  the 
initial  state  | a>  to  all  final  states  | b>  is  : 

W(x,A,  if  1 5)  akafy/jz?  -Z  Fob  (A,^  2^) da9 

b 

*  K  —  ,  /4L/wAm  x0-x)  77?ZpXpJZ^.  pa? 

1(2t)<7?£  *  C  5-1 

where  x=k/k^.  The  squared  matrix  element  ^  2  is  a  lengthy 
express  ion  involving  the  weak  coupling  constants,  V ,  k,  x, X 
,  and  products  of  nuclear  matrix  elements  of  the  form 
described  by  Eg.  4-5.  77£2  is  presented  explicitly  in 
Appendix  C. 

As  enumerated  by  Eood  and  Tolhcek  (Bo  65) ,  the  most 


important  observables  are  : 
1.  the  photon  spectrum 
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5-2 


N(x) dx  =  Z  JJdrt/t  lV(x,A}  2?  $)  dx 


2.  the  photon- neutrino  directional  correlation 

W(x,  8?v)  -  Z  J  w(  x,A,  v,  X,  £) 


3.  the  circular  polarization  of  the  emitted  photons 

Px  m  =  If  da*  da?[w(x,k,  7)-b/(x,A,  r,  X*-/,  s)] 

Jf  da*  da?  [k/(x,A,v,2“+/>  5)  +  h/(x,A,  s)]  s-<* 


4.  the  photon-muon  polarization  directional  correlation 


W(x,  9Ay)  “Z  da?  W(x,d  2,  S) 

X  v 


5-5 


In  practice  the  measurement  of  the  photon-neutrino 
directional  correlation  is  the  most  challenging.  Conceivably 
for  light  nuclei  the  photon  and  the  nuclear  recoil  could  be 
observed  in  coincidence,  although  the  possibility  of  neutron 
emission  from  the  nucleus  would  complicate  matters 
considerably.  At  present  the  photon  spectrum  and  the 
photon-muon  polarization  correlation  have  teen  investigated 
experimentally,  so  that  our  attention  is  focussed  on  the 
computation  of  these  guantities. 
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E.  Photon  Spectrum 

The  differential  photon  spectrum  has  previously  been 
calculated  by  Rood  and  Tolhoek  (Ro  65)  using  the 
closure-har monic  oscillator  model  and  by  Fearing  (Fe  66) 
using  the  giant  dipole  resonance  model-  These  authors 
performed  Pauli  reductions  of  essentially  the  same 
relativistic  Hamiltonian  used  in  this  work,  keeping  terms  in 
the  effective  Hamiltonian  through  0(1/m).  Here  all  terms 
through  0  (1/m2)  arising  from  a  Foldy-Wouthuysen  reduction  of 
the  relativistic  Hamiltonian  of  Eg-  2-1  were  retained, 
resulting  in  the  effective  Hamiltonian  which  appears  in 
Appendix  E.  The  effective  Hamiltonian  was  squared  as 
described  in  Chapter  III,  and  numerical  work  was  done  using 
the  FORTRAN  program  NEW  EMC,  which  evolved  from  a  program 
written  by  H.W.  Fearing.  Cur  photon  spectra  were  calculated 
using  a  squared  matrix  element  complete  through  0{1/m2) •  In 
order  to  facilitate  the  comparison  of  results,  the  following 
set  of  parameters  has  been  defined  for  40Ca  : 
b=2. 03  f m. 

V?  =  85  Mev- 

km= 87- 6  Mev.=  1.03  x  V 

gs=0.0 

gr  =  0.  0 

g  =3.  7 

g^  =  1. 0 

gA  =-1  •  25 

gp=-E.75=  7  x  gA 
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This  set  will  be  referred  to  as  the  standard  parameter 
set  since  it  is  composed  mainly  of  parameters  which  have 
been  used  most  often  in  previous  EMC  calculations 
(Eo  65,  fe  66,  Eo  74).  It  should  be  noted  that  the 
parameters  in  the  standard  set  are  not  necessarily  the  best 
ones.  Choices  for  all  weak  couplings  except  qF  do  however 
reflect  the  current  state  of  experimental  knowledge  of  these 
quantities  (Kr  75).  For  qp ,  the  Goldberger-Treiman  value  g^  =  7xg^ 
was  used.  The  oscillator  parameter  b  was  chosen  to 
reproduce  the  experimental  EMS  charge  radius  for  40Ca 
(Lu  63) .  The  choices  for  the  average  neutrino  energy  V  and 
the  average  maximum  photon  energy  km  were  somewhat  arbitrary 
although  they  tend  to  follow  historical  precedent 
(Ro  65,  Fe  66).  The  ratio  Wt^l-03  is  slightly  larger  than 
the  Rood  and  Tolhoek  estimate  of  1.02  obtained  for  160  by 
equating  results  of  a  closure  calculation  with  those 
obtained  by  summing  partial  transitions,  and  smaller  than 
cur  own  estimate  for  A0Ca  of  1.06,  based  on  somewhat 
different  considerations  explained  in  Chapter  VI.  In 
instances  where  a  relative  rate  (ratio  of  the  rate  for 
radiative  to  ordinary  capture)  is  presented,  the  ordinary 
capture  rate  was  calculated  using  the  effective  Hamiltonian 
cf  Fujii  and  Primakoff  (Fu  59,  ?r  59).  Terms  of  0(1/m2)  in 
the  ordinary  rate  were  not  included  since  these  have  been 
shown  to  contribute  <5%  (Fr  66,  Oh  6  6)  - 

Figure  2  shows  our  0(1/m2)  photon  spectrum  and  a 
corresponding  spectrum  derived  using  the  squared  matrix 
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Contribution  of  previously  neglected  0(1/m2)  terms 
in  the  sguared  matrix  element  to  the  photon  spectrum 
for  4°Ca,  calculated  in  the  CHO  model 
a  -  our  C ( 1/m2)  result 

b  -  using  the  sguared  matrix  element  from  Fearing 
(Fe  66)  - 


lig.  2. 


element  from  Fearing  (Fe  66),  both  of  which  were  calculated 
with  the  closure-harmonic  oscillator  model-  Note  that  for 
k>30  Mev  the  additional  0(1/m2)  terms  contribute  only  a  few 
percent  to  the  spectrum.  This  indicates  good  convergence  of 
the  non- rel ati vistic  expansion  cf  the  Hamiltonian  in  powers 
of  (1/m) ,  as  it  has  also  been  confirmed  that  the  higher 
order  terms  are  individually  small. 

The  contribution  to  the  photon  spectrum  from  various 
parts  of  the  squared  matrix  element  has  been  examined  in 
some  detail  by  Rood  and  collaborators  (Ro  65, Ro  75)  in 
efforts  to  deal  only  with  the  numerically  significant  terms. 
3he  intent  of  the  present  study  has  been  to  retain  all  terms 
where  this  was  possible;  nevertheless  it  remains  instructive 
to  identify  specific  contributions  to  the  photon  spectrum. 
Figure  3  shows  the  photon  spectrum  computed  in  the 
closure-harmonic  oscillator  model  keeping  terms  (including 
velocity  terms)  through  0(1),  G(1/m),  and  0(1/m2),  and  that 
fortion  of  the  rate  through  0(1/m2)  coming  solely  from  the 
velocity  terms,  ie.  those  containing  an  explicit  P.  The  most 
numerically  significant  0(1/m2)  terms  come  from  the  square 
cf  0  ( 1/m)  terms  in  the  effective  Hamiltonian.  These  larger 
0(1/m2)  terms  have  been  included  previously  (Ro  65).  In 
figure  4  the  contribution  of  the  muon  radiating  diagram  and 
its  interference  with  all  other  diagrams  (fig.  1  b-e)  is 
depicted.  As  has  been  known  for  some  time,  the  contribution 
to  the  photon  spectrum  from  the  muon  radiating  diagram  is  a 
major  one  (Pr  59).  Figure  5  shows  the  influence  of  the 
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fig.  3.  Decomposition  of  the  photon  spectrum  for  40Ca 
according  to  powers  of  (1/m)  appearing  in  the 
sguared  matrix  element,  and  the  contribution  of  the 
velocity  terms,  calculated  in  the  CHO  model 
a  -  through  0(1/m2) 
b  -  through  0(1/m) 
c  -  through  0(1) 

d  -  velocity  terms  through  0(1/m2). 
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Fig.  4.  Contribution  of  the  muon  radiating  diagram  and  its 

interference  with  all  other  standard  theory  diagrams 
to  the  photon  spectrum  for  40Ca,  calculated  in  the 
CHO  model 

a  -  complete  0(1/m2)  result 

b  -  muon  radiating  diagram  plus  interferences 
c  -  muon  radiating  diagram  alone. 


10  20  30  40  50  60  70  80  90 

k  (Mev) 


Contribution  to  the  4°Ca  photon  spectrum  from  the 
weak  magnetism  and  induced  pseudoscalar  couplings 
a  -  complete  0(1/m2)  result,  standard  couplings 

b  "  %  =  ° 

c  -  g^=0 
a  -  %=9r  =  ° 


Pig.  5. 
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couplings  and  g^.  Notice  that  the  g^  terms  hecome 
relatively  more  important  as  the  phcton  energy  increases 
while  the  terms  do  not.  The  diagrams  where  the  anomalous 
magnetic  moment  cf  the  proton  and  neutron  represent  the 
source  of  radiation  give  a  net  contribution  depicted  in 
figure  6.  for  these  latter  diagrams  the  contribution  to  the 
rate  is  small  and  they  are  often  neglected.  It  is  estimated 
from  figure  5  however  that  dropping  these  diagrams 
corresponds  to  a  variation  in  qp  of  about  ±  2qA  at  any  given 
photon  energy  k.  From  this  viewpoint  the  small  contributions 
are  seen  to  be  significant  and  hence  should  not  be  dropped. 

The  advent  of  recent  total  photoabsorption  cross 
sections  for  medium  mass  nuclei  (Ah  75)  warranted  fresh 
computations  using  the  giant  dipole  resonance  model*  In  the 
course  of  examining  the  existing  application  of  the  model  to 
radiative  muon  capture  (fe  66),  two  major  improvements  were 
made.  First,  it  was  noticed  that  an  inconsistent  choice  of 
the  average  neutrino  energy  in  the  latter  study  had  resulted 
in  a  giant  dipole  resonance  model  estimate  for  the  relative 
spectrum  which  was  too  low.  Allegedly  the  values  y?  =  85  Mev 
and  3^=68  Mev  had  been  used  in  closure-harmonic  oscillator 
and  giant  dipole  resonance  model  calculations  of  the 
relative  photon  spectrum,  however  the  giant  dipole  resonance 
model  calculation  of  the  ordinary  rate  was  not  in  keeping 
with  the  above  value  for  y.  The  main  nuclear  matrix  element 
sguared  in  the  giant  dipole  resonance  model  appears  as  z 
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.  Effect  of  neglecting  the  anomalous  magnetic  moment 
of  the  proton  and  neutron  on  the  photon  spectrum  for 
40Ca#  calculated  in  the  CHO  model 
a  -  complete  0(1/ffl2)  result 
1  K/o  =  Km  =  0  • 


Fig.  6 


0)0y  =■  / Fit  I  (Mv)uv  ~f~  (Mv)dm 
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Ihe  value  1-67  for  was  taken  directly  from  Foldy  and 

Nalecka  (Fo  64) ,  who  in  fact  used  y=89.8  Mev.  As  a  result  (M 
and  the  ordinary  capture  rate  in  the  giant  dipole 
resonance  model  were  overestimated  as  compared  to  the 
closure-harmonic  oscillator  result,  hence  the  relative 
capture  rate  in  the  giant  dipole  resonance  model  was 
underestimated-  For  v>  =  85  Mev  (M2)^=1.26,  giving  a  value  for 
the  full  matrix  element.  Eg-  5^6,  of  2-28  compared  to  the 
value  2.70  used  in  (Fe  66).  Secondly,  the  new 
photoabsorption  data  for  40Ca  (Ah  75)  permitted  the  removal 
of  an  assumption  made  by  Fearing  concerning  the  relation 
between  the  then  unmeasured  photoproton  cross  section  and 
the  experimental  photoneutron  cross  section  of  Baglin  and 
Spicer  (Ba  64).  It  had  been  assumed  that  the  shapes  of  these 
cross  sections  were  the  same  so  that  the  photoneutron  cross 
section  was  scaled  to  give  the  correct  integrated  cross 
section  for  protcn  photoproduction-  Thus  the  relatively 
small  low  energy  part  of  the  (J',  p)  cross  section  below  the 
ij'-i)  threshold  was  neglected-  Our  calculations  using  the 
Ahrens  data,  however,  show  that  this  part  of  the  cross 
section  is  significant  in  the  giant  dipole  resonance  model 
because  the  low  energy  part  of  the  pho toabsorption  cross 
section  is  the  most  strongly  weighted  in  the  photoabsorption 
integrals  Q,  and  Qz  of  Eg-  4-19. 

In  figures  7  and  8  are  shown  the  results  of  making  the 
two  aforementioned  improvements.  Figure  7  gives  the 
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lig-  7.  Integrals  Q,  and  Qz  over  the  experimental 
photoabsorption  cross  section  for  40Ca 
a  -  Q,  ,  data  of  Ahrens  et  al  (Ah  75) 
t  -  Q,  ,  Baglin  and  Spicer  data  (Ba  64) 
c  -  Qz ,  data  of  Ahrens  et  al  (Ah  75) 
d  -  Qz  ,  Baglin  and  Spicer  data  (Ba  64)  . 
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8,  The  relative  photon  spectrum  for  40Ca,  calculated  in 
the  GDE  model  using  the  squared  matrix  element  from 
Fearing  (Fe  66) -  The  ordinary  capture  rate  A  is  also 
giv  en 

a  -  CHO  model  result ; /f  =4.  36x1  06/s 

b  -  using  the  photcabsorption  data  of  Ahrens  et  al  (Ah  75) 
A  =3.  58x  1 06/s 

c  -  Fearing* s  result  using  the  Baglin  and  Spicer 
photoabsorption  data  (Ba  64)  and  the  average  energies 
V  =89-  8  Mev,  k^=8  8  Mev  ;  A  =2.  98x  10 */s 

d  -  using  the  Baglin  and  Spicer  photoabsorption  data; 

A =3. 44x106/s. 
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integrals  Q/  and  Qz  over  the  scaled  Baglin  and  Spicer 
(Ba  64)  data  and  over  the  Ahrens  (Ah  75)  photoabsorption 
cross  sections  fcr  40Ca.  The  Ahrens  data  is  seen  to  be 
responsible  for  a  substantial  increase  in  the  magnitudes  of  Q/ 
and  Qz  .  The  quantity  (MJ)^  appearing  in  the  expression 
for  the  ordinary  capture  rate  (Eg.  5-6)  also  increased  from 
2.56  (Baglin  and  Spicer  data)  to  3.84  (Ahrens  data).  Our 
relative  photon  spectrum  co  nputed  in  the  giant  dipole 
resonance  model  appears  in  figure  8,  along  with  the  curve 
calculated  by  Fearing  (Fe  66) .  In  order  to  best  illustrate 
the  effects  of  the  improvements,  our  calculation  used  the 
matrix  element  derived  by  Fearing.  Other  multipole  and 
velocity  terms  were  calculated  in  the  closure-harmonic 
oscillator  model.  The  important  difference  which  emerges  is 
that  the  20%  reduction  in  the  relative  rate  found  by  Fearing 
for  the  giant  dipole  model  as  compared  to  the 
closure-harmonic  oscillator  model  becomes  very  small  in  the 
present  calculation  of  these  same  guantitres.  For  60<k<80 
Mev,  approximately  40%  of  this  change  in  the  relative  rate 
is  attributable  to  the  emergence  of  better  photoabsorption 
data,  the  remaining  60%  being  due  to  our  more  consistent 
calculation  of  (M poM  ,  leading  to  a  smaller  result  for  the 
ordinary  capture  rate.  Note  however  that  since  different 
choices  of  km  and  V  can  cause  variations  in  the  relative 
rate  which  are  different  for  different  models,  the  striking 
agreement  between  the  relative  rates  calculated  in  the  CHO 
and  GDE  models  for  the  standard  values  k^=87.6  Mev,  ^=85  Mev 


must  be  regarded  as  something  of  a  coincidence-  Thus 
choosing  ^=84  Mev,  ^=79  Mev  (see  figs.  18  and  19)  one  finds 
that  the  relative  rates  in  the  CEO  and  GDR  models  differ 
substantially  for  k>60  Mev-  The  reason  is  that  while  the 
high  energy  part  of  the  GDR  spectrum  is  fairly  insensitive 
to  the  choice  of  and  y7,  the  CHO  spectrum  is  cut  off  at  k=k 
by  the  closure  approximation. 

The  main  advantage  in  using  the  giant  dipole  resonance 
model  to  describe  radiative  mucn  capture  is  that  the 
fraction  of  the  rate  due  to  dipole  transitions  (about  55% 
for  40Ca)  is  directly  evaluated  using  photoabsorption  data. 
The  importance  of  correctly  evaluating  the  dipole 
contribution  is  illustrated  in  figure  9.  The  dipole  part  in 
the  giant  dipole  resonance  model  amounts  to  only  66%  of  the 
dipole  part  in  the  closure- harmonic  oscillator  model,  so 
that  the  absolute  photon  spectrum  is  substantially  reduced 
in  the  giant  dipole  resonance  theory,  resulting  in  much 
tetter  agreement  with  the  experimentally  observed  absolute 
photon  spectrum  for  40Ca. 

The  squared  matrix  element  given  in  Appendix  C  depends 
on  a  large  number  of  variables.  It  will  prove  instructive  to 
examine  the  sensitivity  of  the  photon  spectrum  to  variations 
in  these  variables,  always  with  an  eye  to  the  feasibility  of 
obtaining  information  on  the  weak  couplings.  Variations  were 
carried  out  with  respect  to  the  standard  parameters. 

Efforts  to  detect  the  presence  of  second  class 
(G-parity  violating)  currents  in  nuclear  beta  decay  have  not 
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Fig.  9.  Dipole  contribution  to  the  photon  spectrum  for  40Ca 
a  -  complete  spectrum,  CHO.  model 
b  -  complete  spectrum,  GDE  model 
c  -  dipole  part,  CHO  model 
d  -  dipole  part,  GDE  model. 
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succeeded  in  demonstrating  that  second  class  currents 
actually  exist  (for  details,  see  reviews  Te  77,  Wu  77). 
Should  they  be  found  however,  the  usual  assignment  of  weak 
coupling  constants  with  g5  =  0. ,  gr=0.,  would  have  to  be 
altered,  figure  10  shows  the  relative  photon  spectrum  for 
values  of  the  scalar  coupling  g5=±1.  In  figure  11  the 
relative  photon  spectrum  is  presented  for  values  of  the 
induced  tensor  coupling  gr=±gAf,  with  gM=3.7,  the  standard 
value  predicted  by  CVC  theory. 

As  first  noted  by  Wolfenstein  (Wo  58)  ,  the  photon 
spectrum  exhibits  strong  dependence  on  the  induced 
pseudoscalar  coupling  qp .  figure  12  indicates  the  degree  of 
sensitivity  of  the  relative  photon  spectrum  to  qp,  and  is 
included  here  for  completeness. 

The  importance  of  the  nuclear  model  in  radiative  muon 
capture  has  been  investigated  by  food  and  Tolhoek  (E o  65) . 
They  established  that  the  relative  photon  spectrum  depends 
cnly  slightly  on  the  nuclear  model  used  to  calculate  it.  The 
basis  for  their  conclusion  involved  comparisons  for  several 
approximations  and  models,  including  a  statistical  model  and 
a  closure-harmonic  oscillator  model.  In  the  case  of  harmonic 
oscillator  wave  functions  which  are  parametrized  by  the 
oscillator  parameter  b,  the  relative  photon  spectrum 
calculated  here  was  found  to  be  insensitive  to  modest 
variations  in  the  value  of  b.  This  is  illustrated  in  figure 
13. 

The  food  and  Tolhoek  study  (Eo  65) ,  while  it  utilized  a 
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Fig-  10-  Effect  of  varying  g5  on  the  relative  photon 

spectrum  for  40Ca,  calculated  in  the  CHO  model-  Ihe 
ordinary  capture  rate  A  is  also  given- 
a  -  gs  =- 1 .  ;  A  =3.  48x  1  06/s 

h  -  gs =  0-  ; A  =4.36x106/s 
c  -  gs=+1 .  ; A  =6-57x10 6/s . 
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Fig.  11.  Effect  of  varying  gr  on  the  relative  photon 

spectrum  for  40Ca,  calculated  in  the  CHO  model 
ordinary  capture  rate  A  is  also  given, 
a  -  gr=+g/>1=  +  3.  7  ;  A  =4.  16x1Q6/s 
b  -  gr  =  0.  ; /l  =4.  36x  1 06/s 

c  -  gr=-gA4=-3.7  ;/j  =4.60x106/s. 
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Fig-  12-  Effect  of  varying  qp  on  the  relative  photon 

spectrum  for  40Ca,  calculated  in  the  CHO  model-  The 
ordinary  capture  rate  is  also  given- 
a  -  gp-  10xg^  =  - 1  2-  5  ;  /  =4.  1 3x1  06/s 
b  -  gp  -lxgA  =-  8.  75  ;  /\  =4.36x1  06/s 

c  -  gF  =  i4xgA  =-5  -  0  ;  /)  =4.  65x  1 06/s 

d  -  gF=0.  0  ;  A  =5.11x106/s. 
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fig.  13-  Effect  of  varying  the  oscillator  parameter  b  on  the 
relative  photon  spectrum  for  40Ca,  calculated  in 
the  CHO  model-  The  ordinary  capture  rate  A  is  also 
gi ven - 

a  -  b= 1 . SO  fm  ; A  =4. 0 1 xl 06/s 
b  -  b=2.  1  5  fm  ;  A  =4. 68x10  Vs. 
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Dumber  of  simple  nuclear  models,  did  not  consider  more 
sophisticated  wave  functions  than  those  of  the  simple 
harmonic  oscillator  shell  model.  In  fact  the  Rood  and 
lolhoek  survey  of  nuclear  models  indicates  that  at  least  for 
doubly-closed  shell  nuclei,  different  elementary  nuclear 
models  predict  relative  capture  rates  which  are  in  good 
agreement  with  one  another.  In  our  work  the  investigation  of 
the  nuclear  wave  function  is  taken  one  step  further  by 
computing  the  photon  spectrum  using  the  Hartree-Fock  wave 
function  for  4°Ca  of  Shao  et  al  (Sh  73).  As  mentioned  in 
Chapter  IV,  the  results  for  the  main  nuclear  matrix  element 
squared  given  in  Egs.  4-6  and  4-7  are  to  be  compared  using 
slightly  different  values  for  b.  Ihis  discrepancy  is  not 
crucial  however  since,  as  has  been  demonstrated  earlier,  the 
relative  photon  spectrum  is  quite  insensitive  to  modest 
variations  of  b. 

As  a  test  of  the  SBL  wave  function,  the  charge  form 
factor  : 


(/+  0.055 az)z 
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was  computed  for  40Ca.  Here  F  (q)  is  the  body  form  factor, 
(q)  is  the  nucleon  charge  form  factor,  and  ~q  is  the 
momentum  transfer.  The  nucleon  charge  form  factor  was 
parametrized  as  in  Lim  (Li  73) .  Results  for  the  SBL  and 
simple  harmonic  oscillator  wave  functions,  along  with  an 
experimental  determination  of  the  form  factor  by  Bellicard 
et  al  (Be  67)  are  presented  in  figure  14.  The  closed  form 
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fig.  14.  Elastic  scattering  charge  form  factor  for  4°Ca 

a  -  CHO  model  using  the  SBL  wave  function,  b=1.7Q 
f  m 

t  -  experimental  determination  of  Bellicard  et  al 
(Be  67) 

c  -  CHO  model  using  simple  wave  function,  b=2.03  fm 
d  -  CHO  model  using  SBL  wave  function,  b=1.90  fm. 


expression  for  the  form  factor  in  the  simple  harmonic 
oscillator  model  is  : 
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with  ^-gb,  b-2.  03  fm,  and  using  the  SBL  wave  function  is  ; 

Z ~  (/-  0.346  7£Z  +  0.  0386  7[*-o.  004/6  y*)Z-ixp{- 

5  -9 

with  b=1.90  fm.  The  region  of  three- mo mentum  transfer  p-n  of 
interest  for  radiative  muon  capture  in  40Ca  is  : 


0  —  j  l  ~  0.56  V-'Tri/ 

*  5-10 


from  figure  14  it  is  seen  that  with  the  conventional  choices 
for  b  both  the  simple  harmonic  oscillator  and  SBL  wave 
functions  predict  a  form  factor  which  is  slightly  below 
experiment.  Figure  15  shows  the  absolute  photon  spectrum 
computed  using  SEL  wave  functions,  and  figure  16,  the 
relative  spectrum.  Since  the  means  of  fixing  b  employed  by 
Shao  et  al  (Sh  73)  involved  considerations  rased  on  160  and 
not  *0Ca,  it  was  decided  in  addition  to  fix  b  for  the  SBL 
wave  function  by  fitting  the  electromagnetic  charge  form 
factor  over  the  range  0<j^/b|<0.6  fm-1.  The  fit  gave  b=1.70 
fm,  and  the  associated  absolute  and  relative  photon  spectra 
also  appear  in  figures  15  and  16.  It  should  be  noted  that 
although  the  absolute  spectrum  is  decreased,  the  relative 
spectrum  is  very  nearly  the  same  as  for  b=1.90  fm.  In  fact 
for  the  CHO  model  all  relative  spectra  calculated  here  using 
a  fixed  set  of  parameters  but  different  wave  functions,  ie. 
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Fig.  15.  Photon  spectrum  for  A0Ca  calculated  using  SB1  wave 
function 

a  -  CHO  model,  b=1.90  fm 
b  ~  CHO  model,  b=1 . 70  fm 
c  -  GDE  model,  b=1.70  fm. 
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Pig.  16.  Relative  photon  spectrum  for  40Ca  calculated  using 
the  SBL  wave  function.  The  ordinary  capture 
rate  A  is  also  given. 

a  -  CHO  model,  fc=1.70  fm  ; A  =4 . 05x 1 06/s 

b  -  CHO  model,  b=1.90  fm  \A- 4.63x106/s 

c  -  GDB  model,  b=1.70  fm  ; A  =3.44 x 1 0 6/s- 
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either  the  simple  harmonic  oscillator  or  SBL  wave  function 
for  several  different  choices  of  b,  are  in  quantitative 
agreement  to  within  5%.  Thus  one  is  able  to  conclude  that  in 
closure  calculations  the  relative  spectrum  is  insensitive  to 
the  choice  of  nuclear  wave  function.  This  conclusion 
supports  the  similar  findings  of  hood  and  Tolhoek  (fio  65). 

The  variation  of  the  relative  photon  spectrum  with  km/V 
is  illustrated  in  figure  17.  Although  km  can  be  roughly 
estimated  from  the  experimental  photon  spectrum,  in  practice 
it  is  difficult  to  fix  km  accurately  enough  to  obtain 
precise  information  on  the  weak  coupling  constants.  This 
point  is  dealt  with  in  more  detail  in  Chapter  VI. 

Having  examined  the  sensitivity  of  the  photon  spectrum 
to  most  of  the  parameters,  we  now  present  our  best 
theoretical  prediction  for  the  photon  spectrum  calculated  in 
the  various  nuclear  models,  assuming  the  Goldberger-Treiman 
value  g/e>  =  7xg/?  for  the  induced  pseudoscalar  coupling.  The 
ether  weak  couplings  have  been  assigned  the  current  values 
in  the  standard  parameter  list.  Our  previous  consideration 
of  the  40Ca  charge  form  factor  has  indicated  that  the  values 
b=2.03  fm  for  the  simple  harmonic  oscillator  wave  function 
and  b=1.70  fm  for  the  SBL  wave  function  are  appropriate.  The 
average  energy  parameters  determined  from  the  sum  rule 
analysis  carried  out  irr  Chapter  VI  are  V=19  Mev,  >^=84  Mev- 
For  the  GDR  model  the  new  photoabsorption  data  of  Ahrens  has 
teen  used.  F'igures  18  and  IS  show  the  relative  photon 
spectrum  in  the  different  models  using  these  best 
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Fig.  17.  Effect  of  varying  the  ratio  km/\7  on  the  relative 
photon  spectrum  for  40Ca.  The  average  neutrino 
energy  was  fixed  at  V  =85  Mev.  The  ordinary  capture 
rate  is /I  =h.  36x  106/s. 
a  -  k m/V;  =  1 . 06 
b  -  km/P=1 . 03 
c  -  km/?=1-00 
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Fig-  18-  *oca  relative  photon  spectrum  for  the  best 
parameter  choices 

y  =19  Mev,  1^=84  Mev,  calculated  in  the  CHC  model. 
The  ordinary  capture  rate  A  is  also 
given-  Atxp  =(2-29±0-06)  x106/s-  (Ha  77) 
a  -  SBL  wave  function; A  =3- 13x106/s 
b  -  simple  wave  function;  A  =3 -37x1 06/s , 
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fig.  19*  40Ca  relative  photon  spectrum  for  the  test 
parameter  choices 

y -19  Mev,  k^,=  84  Mev,  calculated  in  the  GDE  model- 
The  ordinary  capture  rate  A  is  also 
given.  Acxp  =(2.29±0. 06) x106/s.  (Ha  77) 
a  -  SBL  wave  function; A  =3- 02x 10 6/s 
1  -  simple  wave  f  u  rc ti  cn ;  A  =3  -  03  xl  0 6/s . 
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parameters*  Ordinary  capture  rates  are  also  given-  Compared 
fcith  experiment,  all  of  the  rates  come  out  too  high  by  as 
much  as  50%  for  the  CHO  model  ordinary  rate  and  a  factor  of 
almost  three  for  the  GDR  model  radiative  rate  for  k> 60  Mev. 
In  the  test  cases  the  ordinary  rate  is  30%  high  for  the  GDR 
model  and  the  radiative  rate  is  60%  high  for  the  CHC  model 
for  k>60  Mev. 

Recently  evidence  for  the  quenching  of  the  axial  vector 
strength  ty  up  to  25%  in  medium  and  heavy  mass  nuclei  has 
emerged  (Br  78,  Os  79,  To  79).  for  maximal  quenching 
g^=-0-9,  g/>=-6. 3=7x0^ .  The  corresponding  relative  photon 
spectrum  is  presented  in  figure  20  for  the  SBL  wave  function 
and  the  test  parameters  given  atove.  4  While  the  relative 
spectrum  is  changed  only  slightly,  both  the  ordinary  and 
radiative  absolute  rates  are  decreased  some  60%.  Such  a 
large  effect  is  in  the  right  direction  to  achieve  agreement 
with  experimentally  measured  rates. 

The  relative  photon  spectrum  is  especially  sensitive  to 

the  value  of  the  induced  pseudoscalar  coupling  gp,  hence 

potentially  can  be  used  to  obtain  information  about  g^>.  A 

complication  arises  in  that  the  relative  spectrum  is  also 

sensitive  to  the  average  maximum  phcton  energy  km  which 

cannot  be  accurately  fixed  in  an  a  priori  manner.  The  usual 

approach  to  this  difficulty,  which  will  oe  followed  here,  is 

to  take  both  gp  and  km  as  adjustable  parameters  in  fitting 

4  Note  that  the  parameters  V?  =  79  Mev,  km=84  Mev  were 
determined  assuming  no  quenching,  and  that  a  revised  sum 
rule  analysis  incorporating  quenching  might  yield  different 
values  of  V  and  km. 


k  (Mev) 


fig.  20-  4<>Ca  relative  photon  spectrum  for  the  axial  vector 
strength  quenched  25%.  The  SBL  wave  function  and 
best  parameter  choices  y>  =  79  Mev,  km  =  84  Mev  were 
used  with  g„=-0.S,  g/P=-6. 3=7xg,, .  The  ordinary 

capture  rate  A  is  also  given 
a  -  CHO  model; A  =2. 05x1 06/s 
b  -  GDR  model; A  =1 -98x106/s. 
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th€  experimental  photon  spectrum.  The  fixed  parameters  were 
set  to  the  test  values  just  discussed,  figures  21  and  22 
show  two  representative  fits  to  the  relative  photon  spectrum 
for  5 7<k<90  Mev  measured  hy  Hart  et  al  (Ha  77)  ,  using  our 
C(1/m2)  squared  matrix  element.  The  figures  correspond  to 
cases  1  and  2  respectively  in  table  I.  Complete  results  of 
all  fits  done  are  given  in  table  I.  Similar  kinds  of  fits 
were  done  by  Hart  et  al  (Ha  77)  using  the  theory  of  Hood, 
Yano,  and  Yano  (Ho  74)  which  includes  effects  on  the 
intermediate  states  of  the  muon  and  the  nucleon  due  to  the 
nuclear  Coulomb  field.  Their  best  fit,  done  setting 
k;7r/V=1.02,  gA=~  1-24,  gave  ^=86.511.9  Mev,  gp  =-8. 13±2. 00, 
the  fitted  value  for  gp  apparently  agreeing  with  the 
Goldberger-Treiman  prediction  of  gp  =  7xg/?  =-  8.  75.  Absolute 
ordinary  and  radiative  capture  rates  in  the  EYY  theory  are 
however  still  some  80%  higher  than  experiment-  As  the 
results  in  table  I  illustrate,  the  quenching  of affects  the 
relative  photon  spectrum  only  slightly;  however  as  pointed 
out  earlier  the  absolute  rates  can  be  decreased  by  as  much 
as  60%.  Thus  if  gA  is  truly  quenched  on  the  order  of  25%  in 
40Ca,  the  EYY  theory  may  be  able  to  predict  both  absolute 
and  relative  radiative  capture  rates  using  a  value  for 
gp  consistent  with  the  Goldberger-Treiman  estimate. 


N(k)/A(V)x106  (Mev1) 
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Fig.  21.  Two-parameter  fit  to  the  reiative  photon  spectrum 
of  Eart  et  al  (Ha  77)  using  the  CHO  model  with  the 
SHO  wave  function  and  k7n/V  =  1.06.  Results  for  the 
fitted  parameters  are  kOT=89.1±0.9  Mev, 
g^  =  (-2.£±4.8)  g*.  The  ordinary  capture  rate 
is  /I  =5.3x10  6/s . 


N(k)/A(V)x106  (Mev1) 
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Fig-  22.  Iwo-parame ter  fit  to  the  relative  photon  spectrum 
of  Hart  et  al  (Ha  77)  using  the  GDR  model  with  the 
SHO  wave  function  and  km/y='\.06.  Results  for  the 
fitted  parameters  are  1^=87-013.0  Mev, 
g/,  =  (-0-3±2.0)  g The  ordinary  capture  rate 
is/1  =3.9x1  06/s- 


72 


Table  I 

Fits  to  the  relative  photon  spectrum  measured  by  Hart  et  al  (Ha  77) 

Case  Model  Wave  g.  quenched  k  /v  k  (Mev)  gD/g.  10-6A(s-1 ) 
function  A  25%  m  ra  p  A 


1 

CHO 

SHO 

no 

1.06 

89.1 

+ 

0.9 

-2.8 

+ 

4.8 

5.29 

2 

GDR 

SHO 

no 

1 .06 

87.0 

+ 

3.0 

-0.3 

+ 

2.0 

3.89 

3 

CHO 

SBL 

no 

1 .06 

89.1 

+ 

1 .0 

-3.5 

+ 

4.5 

5.01 

4 

GDR 

SBL 

no 

1.06 

86.4 

+ 

3.6 

-0.2 

+ 

1.9 

3.75 

5 

CHO 

SBL 

yes 

1.06 

90.5 

+ 

2.4 

-7.7 

+ 

6.3 

3.63 

6 

GDR 

SBL 

yes 

1.06 

86.4 

+ 

3.2 

-1.4 

+ 

3.2 

2.45 

7 

CHO 

SHO 

no 

1.03 

89.4 

+ 

1.4 

-1.1 

+ 

3.0 

5.59 

8 

GDR 

SHO 

no 

1.03 

86.9 

+ 

3.1 

0.8 

+ 

1.5 

4.05 

9 

CHO 

SHO 

no 

1  .00 

89.4 

+ 

1.6 

0.7 

+ 

2.3 

5.99 

10 

GDR 

SHO 

no 

1 .00 

86.8 

+ 

3.0 

1.8 

+ 

1.2 

4.24 

11  ' 

CHO 

SHO 

no 

1  .06 

89.1 

+ 

0.9 

0.8 

+ 

1  .0 

5.01 

4- 

'without  any  0(l/m2)  terms 


Table  I-  Results  for  km  and  gp  obtained  from  two-parameter 
fits  to  the  relative  photon  spectrum  of  Hart  et  al 
(Ha  77) .  The  other  weak  couplings  were  assigned  the 
standard  values  referred  to  in  the  text;  the  value 
for  the  axial  vector  strength  quenched  25$  was 
g^=-0.9.  Also  given  is  the  ordinary  muon  capture  rate 

A  - 
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C*  Photon  Asymmetry 

If  the  incoming  muon  in  radiative  capture  is  polarized, 
then  there  will  exist  a  measureahle  correlation  between  the 
direction  of  the  emitted  photons  and  the  muon  polarization 
vector.  In  particular  the  directional  distribution  of  the 
photons  can  be  expressed  (Bo  65)  : 

N(x)  =  Jz-  N(x)  (/+  X  5-A) 

5-11 

where  N  (x)  is  the  photon  spectrum  for  unpolarized  muons.  The 
photon  asymmetry  is  depicted  schematically  in  figure  23- 
Fearing  (Fe  75)  has  shown  that  for  the  standard  theory  of 
radiative  muon  capture  (Eo  65)  including  all  weak  couplings 
except  gs  and  all  Feynman  diagrams  shown  in  figure  1,  the 
photon  asymmetry  satisfies  : 

0C=  /  +  0( 

5-12 

where  m  is  the  nucleon  mass.  Thus  to  first  order  the  entire 
correction  to  CCt'l  comes  from  the  0(1/m2)  terms  in  the 
squared  matrix  element.  Cne  can  make  the  following 
observations  : 

1.  the  leading  0(1)  terms  in  the  squared  matrix  element  all 
originate  from  the  muon  radiating  diagram 

2.  as  pointed  out  by  Fearing  (Fe  75),  the  muon  radiating 
diagram  and  its  interference  with  all  other  diagrams 
contribute  only  to  0C  =  1  to  all  orders  in  1/m  and  for  all 
six  couplings. 

Hence  the  first  order  0(1/m2)  correction  to  0C*^  comes 
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N(x)  =  1/4nxN(x)[1+K|S!cos0] 


Fig.  23.  Pictorial  representation  of  the  angle  9  entering 

into  the  expression  for  the  photon  spectrum  in  the 
case  of  polarized  muons. 
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entirely  from  the  square  of  0(1/m)  terms  in  the  effective 
Eami ltonian  originating  from  diagrams  other  than  the  muon 
radiating  one.  Similarly  the  second  order  C(1/m3)  correction 
to  AVI  comes  from  the  product  of  0  (1/m)  and  0(1/m2)  terms  in 
the  effective  Hamiltonian  originating  from  diagrams  other 
than  the  muon  radiating  one.  Thus  one  is  able  to  compute  the 
photon  asymmetry  to  0  (1  / ir^ )  from  a  knowledge  of  the 
effective  Hamiltonian  to  G{1/m2).  The  calculated  asymmetry 
is  shown  in  figure  24.  Note  that  the  0(1/m3)  correction 
amounts  to  at  most  a  3%  drop  in  the  0(1/m2)  asymmetry,  a 
decrease  consistent  with  the  collective  magnitude  of  the 
0(1/m3)  terms,  lor  the  giant  dipole  resonance  model  the 
asymmetry  had  not  in  the  past  been  computed.  This  entirely 
new  result  appears  in  figure  25.  Over  the  experimentally 
observed  region  50<k<90  Mev  the  curve  is  only  a  few  percent 
lower  than  the  harmonic  oscillator  result.  Apparently  the 
parts  of  the  squared  matrix  element  independent  of,  and 
depending  on  3  respectively,  change  in  a  similar  fashion 
when  the  nuclear  matrix  elements  are  changed,  indicating 
that  the  asymmetry  depends  little  cn  nuclear  particulars. 

It  should  be  emphasized  that  the  above  discussion  of 
the  photon  asymmetry  has  been  limited  to  calculations  done 
within  the  framework  of  the  standard  theory.  With  the 
inclusion  of  Coulomb  and  strong  nuclear  potentials,  the 
theorem  proved  by  fearing  (Fe  75)  is  no  longer  generally 
valid  (see  Appendix  F  for  a  detailed  discussion  of  this 
point)  so  that  the  presence  of  additional  interactions  can 
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rig-  24.  Photon  asymmetry  for  40Ca  calculated  in  the  CHO 
model 

a  -  0(1/m2)  result 
1  -  0(1/m3)  result 


OCik) 


10  20  30  40  50  60  70  80  90 

k  (Mev) 


fig.  25.  Photon  asymmetry  for  *°Ca  calculated  in  the  GEE 
model  to  0(1/m2) 
a  -  CHO  model  0(1/m2)  result 
b  -  GDE  model 
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affect  the  asymmetry  in  lower  orders  of  1/m  than  the  second- 
An  indication  that  the  asymmetry  is  just  so  affected  is 
provided  by  the  work  of  hood  et  al  (Ro  74) ,  who  considered 
Coulomb  and  strong  nuclear  potentials  in  radiative  muon 
capture-  Iheir  computed  asymmetry  is  slightly  higher  than 
the  standard  theory  result,  so  that  while  it  appears 
Rearing  *  s  asymmetry  theorem  no  longer  holds,  it  also  seems 
that  the  effects  of  exterior  potentials  on  the  asymmetry  are 
small  for  *°Ca. 

The  photon  asymmetry  is  virtually  unaffected  by  changes 
in  some  of  the  variables  on  which  it  depends.  Among  the 
variables  in  this  category  are  the  oscillator  parameter  of 
the  simple  shell  model  wave  function  or  the  Hartree-Rock 
wave  function  of  SBI  where  it  is  used,  the  average  maximum 
photon  energy  km,  and  the  velocity  terms.  In  summary  the 
asymmetry  does  not  depend  much  on  the  nuclear  physics-  This 
point  is  well  illustrated  by  the  fact  that  our  best 
predictions  for  the  asymmetry,  corresponding  to  the  best 
parameter  choices  made  in  the  previous  section  for  the 
photon  spectrum,  are  essentially  irdisting  ui  shable  unless 
the  distinction  is  made  between  the  CHO  and  GDR  models.  Thus 
our  best  CHO  model  predictions  both  for  the  simple  harmonic 
oscillator  and  §EL  wave  function  differ  from  the  asymmetry 
in  figure  24  by  less  than  and  similarly  for  the  GDR 

model  asymmetry  in  figure  25. 

The  asymmetry  does  depend  to  a  considerable  degree 
however  on  the  values  of  the  weak  coupling  constants.  In 


examining  this  dependence,  the  asymmetry  will  be  calculated 
to  0(1/m2)  in  the  closure-harmonic  oscillator  model  using 
cur  standard  set  of  parameters.  5  Figure  26  shows  the 
asymmetry  for  g„  =  0;  g/7  =  g/P=0;  g^  =  -  0 . 9 ,  g/,=- 6.  3=7xg/?  ; 
illustrating,  in  view  of  the  gp  presented  below, 
that  OC  exhibits  strong  gA  dependence.  The  latter  value  of  gA 
corresponds  to  the  axial  vector  strength  quenched  25$. 
Figure  27  depicts  the  asymmetry  calculated  for  non-zero 
values  of  g5  and  gr,  g5=±1,  gr=±gyM.  In  these  cases  OC  varies 
at  most  ±10%.  Finally  figure  28  shows  the  photon  asymmetry 
for  values  of  the  induced  pseudoscalar  coupling  g/»  =  0,  ^xg^, 
and  10xg^.  With  g/>  =  0,  OC  assumes  a  value  close  to  1,  while  as 
takes  on  increasingly  negative  values,  OC  decreases  fairly 
rapidly.  Thus  experimental  measurements  of  the  photon 
asymmetry  may  prove  useful  in  fixing  gp  ,  provided  the  other 
couplings  can  be  independently  and  accurately  determined.  It 
would  seem  gA  presents  the  biggest  problem  in  this  regard. 
Hart  et  al  (Ha  7  7)  have  measured  OC  f  or  photon  energies 
57<k<90  Mev,  obtaining  an  average  valu  e  0C=+Q .  90±0. 50.  While 
this  result  is  in  agreement  with  the  standard  theory  and 
with  our  results,  a  more  precise  measurement  of  the 
variation  of  OC  with  k  is  clearly  desireable. 


5  As  fig.  24  seems  to  indicate  that  the  0(1/m3)  terms  change 
the  asymmetry  by  a  few  percent  at  most,  these  terms  have 
been  omitted  here  in  the  interest  of  computational 
expediency.  It  has  been  verified  for  the  cases  g^  =  0,  g^=g/»=0 
discussed  below  that  inclusion  of  the  0(1/m3)  terms  indeed 
affects  the  asymmetry  by  no  more  than  3%. 
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Fig.  26.  Effect  on  the  photon  asymmetry  for  40Ca  of  the  weak 
couplings  qv  and  qA ,  calculated  in  the  CKO  model 
a  -  complete  0(1/m2)  result 
b  -  gv  =0 

c  -  qA  =  0 . 

d  -  g^-0.9,  qp=l*qA  (qA  guenched  25%). 
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lig.  27.  Effect  on  the  photon  asymmetry  for  4  °Ca  of  varying  g$ 
and  gr,  calculated  in  the  CHO  model 
a  -  gr=-gM  =  -3.7 
t  -  g$  =- 1-  0 
c  -  gT  =  +g*  =  +3.  7 
d  -  g^  =  + 1  -  0 
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Fig,  28.  Effect  on  the  photon  asymmetry  for  40Ca  of  varying 
,  calculated  in  the  CHO  model 
a  -  g^  =  0. 0 
b  -  g/,  =  4xg>}  =-5. 0 
c  -  gp~  7xgF  =-  8-  75 
d  -  qp  -  10  xgy?  =  -  1  2.  5 


¥1. 


A  MODIFIED  CLOSURE  CALCULATION 


A.  Introduction 

The  calculation  of  radiative  muon  capture  rates  in 
nuclei  has  teen  pursued  by  a  number  of  authors  with  the  aim 
of  extracting  the  magnitude  of  the  weak  induced  pseudoscalar 
coupling  gp  from  measurements  of  the  differential  photon 
spectrum  or  the  photon  asymmetry.  In  most  of  these 
calculations  (e.g.  Pr  59,  Eo  65)  the  closure  approximation 
on  the  final  nuclear  states  has  been  used  to  express  the 
rate  in  terms  of  the  ground  state  expectation  value  of  a 
two-body  operator.  In  this  approach  an  average  maximum 
photon  energy  km  corresponding  to  an  average  excitation 
energy  of  the  residual  nucleus  is  introduced.  The  rate 
depends  very  strongly  on  km  and  hence  uncertainty  in  km 
poses  a  problem  when  one  attempts  to  extract  gP  from 
experiment.  The  situation  is  entirely  analogous  to  that  for 
ordinary  muon  capture,  where  the  result  of  a  closure 
approximation  calculation  is  very  sensitive  to  the  choice  of 
the  average  neutrino  energy  y . 

Eecent  calculations  (Do  72,  Be  73A,  Be  73C,  Go  74, 

Ko  76)  of  ordinary  muon  capture  rates  employing  sum  rules 
have  succeeded  in  largely  eliminating  the  V  dependence  of 
the  rate.  In  particular,  Bernabeu  (Be  73A)  has  shown  that  a 
first  order  expansion  of  the  capture  rate  about  y  gives  a 
result  which  is  essentially  independent  of  the  specific 
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value  of  y  within  a  range  of  plausible  values. 

In  the  present  chapter  these  ideas  are  extended  to 
radiative  muon  capture  and  it  is  shown  that  a  first  order 
expansion  of  the  expression  for  the  photon  spectrum  about  km 
gives  a  corrected  expression  for  the  spectrum  which  is 
icuch  less  dependent  on  the  specific  value  of  km  used  than 
the  usual  closure  result. 

In  Sec.  B  a  corrected  expression  for  the  differential 
photon  spectrum  is  derived.  As  for  ordinary  capture,  two 
additional  terms  arise  besides  the  usual  closure  term.  One 
can  be  calculated  from  the  closure  term  and  the  other  one  is 
evaluated  using  a  generalized  Ihomas-Reiche-Kuhn  { TRK)  sum 
rule  (Er  70) . 

Several  further  corrections  are  considered  which  were 
not  included  in  earlier  calculations  of  ordinary  capture.  In 
particular,  the  effect  of  the  Coulomb  energy  difference 
between  initial  and  final  nuclear  states  is  included  in  a 
simple  way.  Exchange  effects  are  also  included,  albeit  in  a 
phenomenological  way,  by  modifying  the  sum  rule  term  by  an 
overall  factor  which  is  obtained  from  a  sum  rule  analysis  of 
total  photo-absorption  cross  sections. 

In  Sec.  C  the  meaning  of  these  average  quantities  km  and 

V  is  clarified  and  it  is  shown  that  by  a  simple 
extension  of  the  sum  rule  technique,  a  consistency  equation 
which  allows  an  explicit  calculation  of  average  quantities 
related  to  km  and  V  can  be  obtained.  Then  as  an  example  in 
Sec.  D  our  results  are  applied  to  the  closed  shell  nucleus 


using  a  simple  harmonic  oscillator  model,  and  finally 
a  trief  discussion  of  these  results  and  our  conclusions  is 
presented . 
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E.  Theory 

In  an  approximation  which  neglects  the  "velocity 
terms",  the  differential  photon  spectrum  for  radiative  muon 
capture  can  be  written  : 


N(k)  -  &Z/4nz)  /pzZ  I h,  2) 


6-1 


Here  ^  is  the  muon  wave  function,  and  the  sum  on  X  is  over 
the  circular  polarizations  of  the  photon  emitted  with 
absolute  value  of  momentum  k.  The  function  I2(k,A)  contains 
the  nuclear  matrix  elements.  If  one  assumes  for  radiative 
capture  the  relations  given  in  Eg-  3-11,  then  : 


/ 


In  this  expression  y  =  k-V  and  s^  =  (k  +  y  •  with  y^  ,  k^ 
ah 

and  km  respectively  the  neutrino  momentum,  the  photon 
momentum,  and  the  maximum  photon  energy  corresponding  to  the 


transition  a-*b.  G2  is  a  function  of  km  ,  y,  and  the  weak 
coupling  constants,  and  the  sum  on  ab  denotes  an  average 
over  initial  and  a  sum  over  final  nuclear  states.  The 
maximum  energy  available  to  the  photon  is  given  by  : 
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aJ> 


Am  *  777^  '  C 771*- 77l/>)  -Ebb  -  (£7>-£b) 
s  £-(&-£*) 


6-3 


in  which  IB£  is  the  muon  binding  energy  and  Etf ,  E ^  are  the 
energies  of  the  nuclear  states. 

Since  G^2  depends  only  weakly  on  through  the 

combination  k^f /2m  ~  0.05,  we  replace  in  G2  by  an 

&b 

appropriate  average  value  k^,  and  define  : 

/ 


In(k)=J  dyyn7Llai>(Am,k) 


_/  -  -  6-4 

Ihus  the  differential  photon  spectrum  can  be  written  : 


/VC A)  9  (cc77y  G-Z/47c2)  / fy,!* Z,  Cn  In  (A) 

n 


6-5 


The  Cn  ,  which  are  obtained  by  extracting  the  explicit  y 

dependence  from  G2  ,are  functions  of  k—  and  the  weak 

ah  771 

couplings.  When  the  In  are  evaluated  in  the  closure-harmonic 

oscillator  model.  Eg.  6-5  gives  an  expression  for  the 

spectrum  in  closure  approximation,  N{kw,k),  which  exhibits 

strong  dependence  on  the  average  maximum  photon  energy  km, 

as  will  be  seen  in  the  next  section. 

In  an  attempt  to  remedy  this  situation,  correction 

terms  to  the  closure  approximation  were  calculated  using 

techniques  analogous  to  those  used  by  Bernabeu  (Be  73A)  for 

oJ> 

ordinary  capture.  Thus  1^  ( km  ,k)  was  expanded  to  first 
order  in  k°£  about  an  average  value  km  and  a  corrected 

t»  <xb 

expression  I2^  for  which  the  km  dependence  is  explicit  and 
linear  was  obtained. 

I tii  =  Iai  (km,  k)J^^m  +  CkUt  (k™^k)//£°k».  (km  ~  km) 


87 


Using  Eg.  6-3  : 

IaM  ~  (/+  (£~ km)  3km  )ZdJ>  (km.  A)  ~  9km  £(Zi>  ~£n)I (km^A)]/ ^ 

6  -6 

Substituting  Eg.  6-6  in  Eg.  6-4,  a  new  expression  is 
obtained  : 

/ 

In  fa.  A)  x  (/+(£-kJ§km)ln (km,k) -Jt(yyn?kn  \ZkOu-k)* 

~l  TTJu. 

x  (Ek-Ea.)  8 (kn-k) /<Z/Z  7-  (i)  f  ^  '*/*>/ 2 km  I 

i  J 

6-7 

where  In  (k^k)  is  the  closure  approximation  to  In{k).  When 

substituted  in  Eg.  6-5  this  yields  an  expression  for  the 

differential  photon  spectrum  which  will  be  denoted  N(km,k) . 

The  final  term  in  Eg.  6-7  can  be  calculated  using  a 

modified  Thomas-Beiche-Kuhn  sum  rule  { Er  7  0)  .  For  a 

many-particle  Hamiltonian  H  =  I+V,  an  operator  71 0^  which 

JL 

satisfies  the  commutation  relation  Vj  -  0  then  also 

y 

satisfies  : 

>2  CEn-Ex)  /<A/Z &/*>/*  m  M  <&■! JL(ViOi)’ (V±  OE) I E> 

i  i  Zm  i  6_8 

with  E* ,  E*  eigenstates  of  H  and  m  the  nucleon  mass. 

Ah 

There  are  two  refinements  which  can  now  be  made  before 
applying  this  sum  rule  to  the  evaluation  of  Eg.  6-7.  These 
have  not  generally  been  made  in  previous  calculations  for 
ordinary  capture,  but  do  seem  to  have  numerical  significance 
and  so  will  he  included  here.  The  first  deals  with  the 
Coulomb  energy  shift.  Observe  that  Eg.  6-8  contains  EJ-E^, 
the  difference  in  eigenvalues  of  the  nuclear  Hamiltonian, 


88 


whereas  Eg.  6-6  involves  E^-E^,  the  difference  in  the  actual 
nuclear  level  energies,  which  cf  course  includes  the  Coulomb 
energy.  The  two  are  related  by  E^-E^  =  E^-E^-Ec  where  Ec  is 
the  Coulomb  energy  difference  between  the  (A,Z)  and  (A,Z-1) 
ground  states.  Substitution  of  this  relation  into  Eg.  6-7 
effectively  replaces  E  by  E+Ec  in  the  second  term. 

The  second  refinement  has  to  do  with  the  influence  of 

exchange  corrections.  When  the  TEK  sum  rule  is  applied  to 

photoabsorption  processes  it  fails,  predicting  total 

photoabsorption  cross  sections  which  are  too  low.  This  is 

not  entirely  unexpected  since  it  has  been  known  for  some 

time  that  nuclear  exchange  potentials  not  satisfying 

[ 7L  O+j  Vexch ]  =  0  give  rise  to  terms  which  make  a 

i 

significant  contribution  to  the  sum  rule  (Le  50,  Ee  57).  The 
inclusion  of  exchange  forces  would  enhance  the  right  hand 
side  of  Eg.  6-8  by  an  amount  conventionally  described  by  a 
phenomenological  factor  (1  ♦  £),  which  for  specific  nuclei 
may  be  determined  from  fits  to  experimental  total 
photoabsorption  cross  sections.  Moreover  since  the  processes 
cf  photcabsorption  and  muon  capture  are  both  dominated  by 
dipole  transitions  and  in  that  approximation  involve  similar 
operators,  we  may  reasonably  expect  that  values 
for  S  determined  from  phctoabsorption  data  (Ah  75)  can  be 
used  to  estimate  the  s  necessary  to  calculate  muon  capture 
rates.  Thus  use  will  be  made  of  the  modified  TEK  sum  rule  in 


the  form  : 
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Z  (£t-£L)/<61Z  CL  ln>lz=  &£&<*/£  (Vi  o/)-  (ioi)M> 

b  b  Z.771  £ 

Applying  Eg.  6-9  to  the  final  term  in  Eg.  6-7,  the 
differential  photon  spectrum  obtained  is  : 

N  ( km,  k)  =(cc 77^  GZ/4tTZ)J^!ZZ  Cnln(km,k) 


6-9 


6-10 


with 


In  (km,k)~  [/+  (Et  £c  ~  km)  36m] In  (km,  k)  ~  (!+  S)  AM± 


771  77^ 


X 


f(Zkm.  ~&kkm.  +  7k  km- 3k3)/ (n+j)  71  eYC7l 

(3kk?n  ~ 6k*kj?t  +  3 k3) / (n+z)  n  odd 


6-11 


To  obtain  the  relative  rate,  the  ordinary  capture  rate 
calculated  in  the  same  approximation  is  required.  Thus  the 
ordinary  rate  A  is  written  (again  neglecting  velocity  terms) 


as 


A 


71 


/T/I  7L  (&f  +  3  A/!a6  (Vat>) 


6-12 


with  : 


m*  (v*i)  -  (**/*bf J %  /&izr-(fi£iU'A' '/*>! 

J 


6-13 


Proceeding  as  for  radiative  capture  the  corrected  rate  is 
found  to  be  : 

,Z  ,  2  _  2  \  rsf  z 


AC y)  *  77U  /%*/  (&F  +  3 &&r)  M  c?) 


Zrr 


6-14 


where  : 


rsfZ 

M  (?) 


-  C/t  (£+Et-?)&] M%)-Z(lti)ZV3^(XA 


m 


6-15 
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with  M2(y>)  the  usual  closure  result.  This  is  essentially  the 
result  of  Bernabeu  (Be  73A) ,  with  Coulomb  and  exchange 
corrections  added  as  above. 

So  far  in  this  and  previous  work  the  velocity  terms 
have  been  neglected.  This  has  been  necessary  for  the  sum 
rule  piece  because  an  operator  does  not  satisfy  the 
assumptions  necessary  to  obtain  the  simple  sum  rule  of  Eg. 

Observe  however  that  both  the  velocity  terms  and  the 
sum  rule  corrections  to  the  main  terms  are  0(1/m).  Thus 
presumedly  a  sum  rule  correction  to  the  velocity  terms, 
which  would  involve  commutators  with  the  kinetic  energy 
j?2/2m#  would  be  of  0(1/m2),  i.e.  of  the  same  order  as  other 
terms  neglected.  Corrections  of  0(1/m)  are  obtained  however 
by  including  the  velocity  terms  in  the  usual  closure  results 
7L  Cw  ln  (kw/k)  and  (g2  +  3G2)  M2  {/)  appearing  in  the 

71  •" 

first  two  terms  of  Egs.  6-11  and  6-15.  So  henceforth  the 
velocity  terms  will  be  incorporated  in  £  6-7J,  In  (km#k)  and  (g^2  +  3G^ 
X  M2  (y?)  ,  yielding  a  result  which  is  expected  to 
be  accurate  to  0  (1/m)  . 

C.  Sum  fiule  Gelations  for  Average  Excitation  Energies 

Before  actually  evaluating  the  expressions  derived 

above  in  the  context  of  a  specific  model,  a  discussion  of 

the  meaning  of  the  parameters  k^and  V  is  presented.  In  the 

process,  a  consistency  relation  is  derived  by  using  the  same 

sum  rule  technigues  employed  above.  The  relation  allows  us 

to  make  an  explicit  estimate  of  the  "physical"  values  of  km  and 

V, 


So  far  k^  and  ?  have  teen  somewhat  loosely  referred  to 
as  "average"  maximum  photon  energy  and  "average"  neutrino 
energy  corresponding  to  an  average  nuclear  excitation  of  the 
residual  nucleus.  Strictly  speaking  they  are  not  averages  in 
the  physical  sense  hut  simply  parameters  defined  formally  by 
N(k)  =  N(k7n#k)  and  A  ^  A  (?)  where  N  (km,k)  and  A  (?)  are  the 
closure  approximations  to  the  actual  rates  N(k)  and  A.  Thus  k 
and  V  are  the  values  which,  when  used  to  evaluate  the 
rates  in  closure  approximation,  give  the  correct  answer. 

In  the  previous  section  improved  approximations  to  the 
correct  rates  were  obtained,  i.e.  N  (kw,k)  and  A  (?)  •  It  will 
be  shown  in  the  following  discussion  that  at  least  in  a 
simple  model  these  are  relatively  independent  of  km  and  ? . 
Ihus  the  appropriate  values  of  km  and  V  to  use  in  a  closure 
calculation  can  be  estimated  by  solving  the  equations 

/v 

N(km,k)  =  N(k;w,k)  and  A  (?)  -  A  (?)  *  that  is,  by  determining 
the  intersection  point  of  the  closure  and  corrected 
calculations.  Ihis  is  just  the  point  where  the  two 
correction  terms  in  Eg.  6-11  or  those  in  Eg.  6-15  exactly 
cancel.  The  resulting  parameters  are  then  those  which  must 
be  used  in  a  closure  approximation  calculation  for 
consistent  results.  They  are  of  course  somewhat  dependent  on 
the  model  used  for  the  nuclear  matrix  elements  and  thus  if 
drastically  different  values  are  required  to  fit  the  data, 
one  should  view  the  model  with  suspicion. 

The  "physical"  values  of  average  maximum  photon  energy 
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and  average  neutrino  energy  (which  we  shall  write  as  kw  and 
V  )  are  in  principle  different  from  the  parameters  k^  and 


V.  Their  p^  moments  can  be  defined  formally  as  : 


km  =  Z  (kmf  Nafk)  /Z  /Va  (A) 


Z  (Vaf Aa  /Z  Aa 


6-16 


where  the  denominators  are  just  the  rates  N(k)  and  A,  and 
the  numerators  are  the  appropriate  quantities  for  the 
transition  a  — >  b  weighted  by  the  rates  for  that  transition 
and  summed  over  all  states  b. 

In  closure  approximation  km  =  km  and  V  =  V  ,  which  is 
what  is  normally  assumed.  Corrections  to  these  relations  can 
however  be  calculated  in  exactly  the  same  way  that 
corrections  to  the  closure  approximation  for  the  rates  were 
calculated,  so  that  when  such  corrections  are  included, 
these  equalities  will  no  longer  hold.  To  do  this  let  us 


expand  both  numerator  and  denominator  of  Eg.  6-16  about  km  =  km 
and  Va&  -  V  and  carry  through  the  sum  rule 
evaluation  as  done  in  Egs.  6-6  through  6-11.  The  results, 
keeping  only  the  first  order  correction,  are  : 


krS=  km  A />  km  ' '  f(£+Et' km)  ~(l+l)  (ZA/mvpA  (Arn^f 


7L  Cn  Zn  (km}  k) 


6-17 


and  : 


. 
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Yr=  Vr  +  P9P~'  {(Ei£t-Y)-(l+i)VtZ/(ZmMZM)}  6-18 

Thus  for  p=  1  these  relations  express  the  physical  averages  k m 
and  V  as  the  closure  parameters  km  and  V  plus  a 
correction  term.  Note  that  the  correction  is  essentially  the 
same,  only  without  the  derivatives,  as  that  appearing  in  the 
equations  for  N  (kw,  k)  and  /^(y?)  - 

A  more  useful  result  can  he  obtained  by  expanding  the 
right  sides  of  Eg.  6-16  about  km  and  V  instead  of  km  and  V>. 

The  resulting  equations  (identical  to  those  above  with  km 

and  everywhere)  now  provide  consistency 

relations  which  can  be  solved  for  km  and  ?.  Since  the 
leading  terms  cancel,  the  solutions  correspond  to  the  zero 
cf  the  term  in  brackets  and  provide  an  estimate  of  the 
physical  averages  corresponding  to  the  physical  average 
excitation  energy  of  the  residual  nucleus.  Note  that  the 
results  will  depend  on  the  model  used  for  the  nuclear  matrix 
elements  in  In  (k^,k)  and  M2  [V) •  These  consistency  relations 
will  be  evaluated  for  a  simple  model  in  the  next  section. 

£.  Application  to  *°Ca 

To  illustrate  the  ideas  outlined  in  the  previous 
sections,  an  application  of  our  method  will  be  made  to  the 
nucleus  *0Ca,  the  nuclear  matrix  elements  being  evaluated 
using  harmonic  oscillator  shell  model  wave  functions.  Such  a 
model  is  perhaps  too  simple,  but  has  been  used  for  most 
comparisons  with  data  and  in  any  case  will  illustrate  the 
basic  results.  In  this  model  the  closure  result  for  ln  is  : 


6-19 


/ 


where  for  40Ca  with  ^72  -  (sb)2  the  main  unclear  matrix 
element  is  : 


rf-io  [/-(/+£-  n?+ z£) exp (-n/z)] 

°°  *>40  6-20 

The  oscillator  parameter  is  taken  to  be  b  =  2.03  fm. 
Expressions  6-5  and  6-10  for  the  differential  photon 
spectrum  were  evaluated  usirg  the  Cn  from  the  appendix  in 
fearing  (Fe  66)  and  the  set  of  weak  couplings  gv  =  1.0# 

=  3-7,  gA  -  -1.25,  gp  =  7  gA ,  g^  =  0,  gT  =  0.  From  Nolen 
and  Schiffer  (No  69)  Ec  =  7.13  Mev,  and  from  Engfer  et  al 
(En  74)  EB£  =  1.066  Mev,  which  give  the  value  E  +  Ec  =  110.4 
Mev.  As  noted  above,  the  velocity  terms  have  been  included 
as  has  the  phenomenological  correction  for  exchange  effects, 

+  S)  • 


It  has  been  customary  to  present  results  for  radiative 
muon  capture  rates  as  a  ratio  of  the  differential  photon 
spectrum  to  the  ordinary  rate,  as  presumedly  some  of  the 
model  dependence  of  the  nuclear  matrix  elements  will  then 
cancel,  though  factors  which  independently  affect  the 
overall  scale  of  the  amplitudes  for  radiative  and  for 
ordinary  capture  will  of  course  also  affect  the  ratio.  In 
figures  29  and  30  are  shown  such  plots  which  illustrate  the 
main  features  of  our  results.  It  can  be  seen  that  the  ratio 
of  corrected  quantities  N  (kw,  k)  //{[  y?)  ,  shown  as  solid  and 
short-dashed  curves  for  two  different  values  of  S,  is 


95 


76  80  84  88  92  96  100 


kjMeV) 


fig.  29.  The  relative  capture  rate  N(k)//t  for  4°Ca  for  k=60 
Mev  and  k„,/>>=1.06.  The  usual  closure  result 
(long-dashed  curve)  is  compared  to  the  corrected 
result  f or S  =0  (short-dashed  curve)  and  for  6  =1-15 
(solid  curve)  -  The  horizontal  lines  are 
experimental  bounds  on  N(k)//1  from  (Ha  77). 


(MeV-1) 
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Fig.  30.  The  relative  capture  rate  N(k)//f  for  4°ca  for  k=75 
Mev  and  k7n/);  =  1-06.  The  curves  have  the  same  meaning 
as  in  figure  29. 


97 


generally  much  less  dependent,  on  k^  than  the  usual  closure 
result  N  (k^k )/A{V)t  shown  as  a  long-dashed  curve-  Thus  the 
important  gualitative  result  is  that  for  radiative  capture, 
just  as  for  ordinary  capture  (Be  73A)  ,  the  sum  rule 
technique  allows  one  to  ohtain  a  result  more  or  less 
independent  of  the  closure  parameter  km  over  a  reasonable 
range  of  km.  Note  that  since  A  (V)  is  itself  nearly 
independent  of  V  in  this  approach,  these  qualitative 
features  hold  for  the  absolute  rate  N(km,k)  just  as  for  the 
ratio  N  (kmt  k )  //T ( V7)  - 

It  is  also  observed  that  the  curves  for  the  corrected 
ratio  N  (km,k) /A  (V7)  and  the  closure  ratio  N  (km, k) /A{V) 

fSJ 

intersect  in  a  region  where  N (km,k) /A (/) ,  and  due  to  the 
stability  of  A  (Y)  ,  N(km,k)  itself,  is  stable.  Thus  the 

A/ 

intersection  of  N(km,k)  and  Nfk^k)  can  be  made  a  criterion 
in  the  selection  of  a  value  of  km  to  be  used  in  an  ordinary 
closure  calculation  of  N (k) .  Similarly  the  appropriate  value  of 

r*J 

V  can  be  determined  from  the  intersection  of  A  (v7) 
and  A  (V7)  -  As  emphasized  in  Sec.  C,  km  and  V  determined  this 
way  are  basically  just  parameters  which  force  the  closure 
approximation  to  give  results  for  N  (kw,k)  and  A  (V)  equal  to 
the  sum  rule  corrected  values,  which  in  turn  approximate  the 
correct  results. 

In  figures  29  to  31  are  shown  results  for  two  different 
values  of  £,  5=0  and  £  =  1.15-  The  value  6=0  corresponds  to  no 
exchange  contribution  correction  to  the  sum  rule  piece,  whereas 
5=1.15  is  the  lower  limit  on  6  given  by  Ahrens  et  al 
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fig.  31.  Influence  of  the  Coulomb  energy  shift  and  velocity 
terms  on  the  40Ca  relative  capture  rate  for  k=70 
Mev  and  km/>>=1.06.  Eesults  including  velocity  terms 
(solid  curves)  and  not  including  them  (dashed 
curves)  are  shown  for  two  choices  of  the 
phenomenological  factor  6  and  with  and  without  the 
Coulomb  energy  shift  Ec .  The  horizontal  lines  are 
experimental  bounds  on  N  (k)//J  from  (Ha  77) 
a  -  usual  closure  result 
b  -  5  =  1.15,  Ec  =7.  1  Mev 
c  -  5=1.15,  Ec  =0 

d  -  5  =  0 ,  Ec  =7.  1  Mev 
e  -  5  =  0,  Ec  =0. 
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(Ah  75)  as  determined  from  the  ratio  of  the  experimental 
photoabsorption  cross  section  to  the  classical  dipole  sum. 

It  is  observed  that  the  introduction  of  the  correction  &  is 
numerically  important,  particularly  for  large  k.  The 
correction  is  also  important  in  principle  since  a  consistent 
treatment  of  both  photoabsorption  and  radiative  muon  capture 
using  sum  rules  requires  that  the  information  available  from 
photoab sorption  be  used  to  constrain  the  parameters  for 
radiative  muon  capture.  It  seems  reasonable  to  take  the 
values  for  £,  1.15  <  &  <1.4,  directly  from  Ahrens  et  al 
(Ah  75)  ,  since  the  operators  in  the  nuclear  matrix  elements 
for  photoabsorption  are  essentially  the  same  as  for  muon 
capture  (Be  73C,Fo  64)  and  since  for  4 °Ca  the  dipole 
transition  dominates  both  ordinary  and  radiative  capture. 

For  the  minimum  correction  S  =1.15  the  absolute  rates  are 
in  much  better  agreement  with  the  data  (Ha  77)  than  for  no 
correction  (cf.  Table  II).  As  6  is  increased  from  1.  15  to 
1.4,  the  relative  rate  decreases,  in  most  cases  improving 
the  agreement  with  the  data.  The  current  experimental 
uncertainty  in  S  corresponds  to  an  uncertainty  in  km, 
determined  from  the  intersection  criterion,  of  about  1.5 
Hev. 

Figure  31  shows  in  addition  the  sensitivity  of  the 
results  to  the  two  other  corrections  which  have  been 
included.  Coulomb  energy  shifts  and  velocity  terms.  The 
Coulomb  correction  is  very  important.  Its  inclusion  tends  to 
increase  N  (k7n,k ) /A  (V*)  and  the  radiative  rate  N  (km,k)  by  as 
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Table  II.  Ordinary  and  radiative  mucn  capture  rates  A  and  N, 
and  correspon din g  average  parameters  V  and  km, 
obtained  from  the  intersection  of  the  closure  and 
corrected  rate  expressions  for  4  °Ca  in  a  simple 
harmonic  oscillator  model.  Also  presented  are 
values  of  the  physical  averages  V  and  k**,  obtained 
by  solving  the  consistency  relations  (Egs.  6-17 
and  6-18)  for  these  quantities.  The  cases  are  : 

I:  £=1.15,  Efc  =  7 .  1  Mev ,  velocity  terms  included; 

II:  6=1.15,  Ec  =  7 .  1  Mev,  no  velocity  terms; 

III:  6=1.15,  Ec=0-0  Mev,  velocity  terms  included; 

IV:  6  =  0.0,  Ec=7.1  Mev,  velocity  terms  included; 

V:  6=0.0,  Ec=7.1  Mev,  no  velocity  terms. 
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Table  II 

Ordinary  and  Radiative  Muon  Capture  Rates  Obtained  from 

the  Intersection  Criterion 


k 

(MEV) 

V 

(ME  V) 

V 

(MEV) 

km 

(MEV) 

km 

(MEV) 

km^v 

km/v 

10"  7  A 
(SEC"1) 

N 

(SEC-1 

MEV"1) 

106N/A 

(MEV-1) 

1 

60 

78.72 

82.76 

8*1.08 

84.93 

1.068 

1.026 

0.333 

6.55 

1.97 

65 

83.92 

84.42 

1.066 

1.020 

4.96 

1.49 

70 

83.53 

83.73 

1.061 

1  .012 

3.15 

0.95 

75 

82.75 

82.81 

1.051 

1.000 

1.32 

0.40 

1 1 

60 

77.12 

80.78 

82.22 

83.30 

1 .066 

1.031 

0.283 

5.24 

1.85 

65 

82.31 

83.02 

1.067 

1 .028 

3.96 

1.40 

70 

82.17 

82.57 

1.065 

1.022 

2.48 

0.88 

75 

81.74 

81.89 

1  .060 

1  .014 

0.98 

0.35 

1 1 1 

60 

73. 'll 

76.94 

77.86 

78.22 

1.061 

1  .017 

0.259 

3.65 

1.40 

65 

77. 44 

77.54 

1  .055 

1 .008 

2.24 

0.86 

70 

76.67 

76.64 

1 .044 

0.996 

0.83 

0.32 

75 

75.80 

75.51 

1.033 

0.981 

0.01 

0.00 

IV 

60 

93.12 

96.05 

96.99 

97.94 

1.042 

1.020 

0.599 

15.98 

2.67 

65 

97.04 

97.81 

1  .042 

1.018 

13.96 

2.33 

70 

97.04 

97.62 

1.042 

1.016 

11.73 

1.96 

75 

96.99 

97.38 

1  .042 

1.014 

1  t 

9.23 

1.54 

V 

60 

92.31 

95.12 

95.91 

97.00 

1.039 

1.020 

0.540 

13.87 

2.57 

65 

96.01 

96.96 

1  .040 

1.019 

12.20 

2.26 

70 

96.  14 

96.87 

1  .041 

1.018 

10.33 

1.91 

75 

96.18 

96.71 

1.042 

1  .017 

8.17 

1.51 
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much  as  a  factor  of  two  or  three  for  some  k  in  extreme 

r*t 

oases.  It  also  increases  A  (V) ,  though  by  a  smaller  amount, 
and  increases  k^,  and  obtained  from  the  intersection 
criterion  by  roughly  the  Coulomb  energy  shift.  Thus  in  this 
approach  it  is  clearly  important  to  consider  the  Coulomb 
energy. 

The  velocity  terms  are  somewhat  less  important,  but 
still  should  be  included.  They  increase  the  ordinary  rate  by 
15-20%,  the  radiative  rate  by  30-35%,  and  the  ratio  by 
5-15%.  This  is  a  somewhat  larger  effect  than  is  found  in  the 
usual  closure  approximation  (where  the  velocity  terms  make 
very  little  difference  in  the  ratio  N  (km,k)  /A{V)  )  both 
because  the  corrected  rates  N(kw,k)  and  A  (V)  are  more 
sensitive  to  the  velocity  terms  than  their  closure 
counterparts  and  because  these  changes  in  the  rates  change 
the  values  of  km  and  V  obtained  from  the  intersection 
criterion  by  about  1  Mev  (cf.  Table  II)  which  in  turn 
produces  a  change  in  the  contribution  of  the  leading  terms. 

finally  the  results  obtained  for  km  and  v7  from  the 
intersection  criterion  and  those  for  the  physical  averages  k^ 
and  V  obtained  by  solving  the  consistency  relations  of 
Zgs.  6-17  and  6-18  merit  some  discussion.  These  results  are 
given  in  Table  II  where  V,  km,  km/V,  the  ordinary  and 
radiative  rates  and  their  ratio  evaluated  at  km  and  V,  and 
the  values  v>#  km  and  km/V  are  shown  for  several  combinations 
of  the  corrections  discussed  above  and  for  several  values  of 
k. 
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The  most  striking  result  is  that  for  fixed  &,  the  ratio 
km/y  is  essentially  constant  even  though  km  and  V  vary 
individually  by  sizeable  amounts  and  for  some  k  the  rates 
^an  vary  by  as  much  as  factors  of  two.  Thus  it  seems 
sensible  to  fit  data  by  fixing  km/¥  and  treating  say  km  as  a 
parameter,  as  was  done  by  Hart  et  al  (Ha  77).  For  the 
standard  Hood  and  Tolhoek  model  (Ec  65)  ,  using  simple 
harmonic  oscillator  wave  functions  and  £=1.15,  it  is  seen 
from  Table  II  that  km/V  =  1.06  ±0.01  is  an  appropriate 
value-  For  the  less  desireable  choice  6=0,  km/?  is  still 
fairly  constant,  but  a  bit  lower,  about  1.04- 

These  values  of  km/y  are  somewhat  larger  than  those 
used  before.  This  can  be  understood  by  looking  at  the  ratio 
of  "physical'1  averages  km/V  obtained  by  solving  the 
consistency  relations  starting  from  the  definitions  of  Eg. 
6-16.  The  values  are  not  quite  as  constant  as  was  the  case 
for  km/V ,  but  still  vary  by  only  ±  00.02.  They  are  generally 
lower  than  km/V .  Values  of  km/V  used  in  previous  fits  to 
data  have  usually  been  obtained  from  calculations  such  as 
that  of  Eood  and  Tolhoek  (Eo  65),  who  got  1.02  for  16C  by 
comparing  the  closure  result  with  that  obtained  by  summing 
over  partial  transitions.  Such  a  calculation  starts  with  the 
definitions  of  "physical"  averages  {Eg.  6-16)  however,  and 
so  is  really  a  calculation  of  km/V~  Hence  the  agreement  with 
our  consistency  result  for  km/v 7  is  encouraging.  The  effect 
of  the  additional  sum  rule  correction  terms  is  to  increase 
slightly  the  value  of  km/V  appropriate  for  a  closure 
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calculation. 

It  should  be  emphasized  that  the  specific  numbers,  e.g. 
for  Ith/V ,  may  depend  on  the  model  used,  in  this  case  the 
harmonic  oscillator  model,  though  the  gualitative  features 
it  is  presumed  are  general.  In  a  mere  complete  calculation  a 
number  of  additional  corrections  should  be  included,  in 
particular  the  propagator  corrections  of  Hood,  Yano,  and 
Yano  (Ho  74) ,  the  higher  order  terms  suggested  by  Fearing 
(Fe  75)  ,  perhaps  better  wave  functions,  and  perhaps  the 
requirement  of  consistency  with  electromagnetic  matrix 
elements  as  done  in  the  GDB  model  (Fe  66,  Fo  64). 

It  is  interesting  to  note  that  with  £=1. 15  the  ordinary 
rate  3.  33  x  106  sec*1  is  in  fair  agreement  with  the 
experimental  values  (2.29  ±  0.06)  x  106  sec-1  (Ha  77)  and 
(2.53  ±  0.02)  x  106  sec*1  (Di  71).  The  radiative  rate  is 
somewhat  high,  but  presumedly  will  be  reduced  on  the  order 
of  20%  by  the  EYY  correction  (Ho  74),  so  that  it  may  be 
possible  to  achieve  gualitative  agreement  here  as  well. 
Alternatively,  an  increase  in  S  (£=1.15  is  a  lower  limit) 
reduces  ordinary,  radiative,  and  relative  rates. 

Thus  to  summarize,  it  has  been  shown  that  the 
differential  photon  spectrum  for  radiative  muon  capture  can 
be  calculated  in  a  way  which  renders  it  nearly  independent 
of  the  value  of  the  maximum  photon  energy,  in  a  fashion 
analogous  to  that  for  ordinary  capture.  The  procedure 
consists  of  expanding  the  expression  for  the  photon  spectrum 
to  first  order  about  An,  then  using  closure  to  perform  the 


sum  over  final  states.  Evaluation  of  one  of  the  resulting 
correction  terms  is  done  using  a  modified  TRK  sum  rule, 
leading  to  an  expression  for  the  photon  spectrum  which 
exhibits  only  slight  dependence  cn  k^.  The  intersection  of 
this  result  with  the  usual  closure  result  determines  values 
of  km  and  for  use  in  closure  approximation  calculations. 
From  a  practical  point  of  view  the  most  useful  result  may  be 
that  in  the  standard  Rood  and  Tolhoek  model  (Bo  65)  the  sum 
rule  correction  gives  values  for  km/V  which  are  essentially 
constant,  and  slightly  larger  than  those  used  previously. 


VII. 


SOJifJABY  AND  CONCLUSIONS 

In  the  present  work  the  *°Ca  photon  spectrum  has  been 
calculated  to  0(1/m2)  and  the  photcn  asymmetry  calculated  to 
C(1/m3)  ,  using  the  impulse  approximation  with  the  standard 
theory  set  of  Feynman  diagrams  shown  in  fig.  1.  While  some 
cf  the  0(1/m2)  terms  previously  considered  (Bo  65)  are  not 
negligible,  the  collective  magnitude  of  previously  neglected 
C(1/m2)  and  0(1/ra3)  terms  was  found  to  be  small,  resulting 
in  a  change  over  the  experimentally  accessible  range  k>55 
Mev  of  less  than  5%  for  both  the  photon  spectrum  and 
asymmetry  as  compared  with  the  standard  theory.  Introduction 
of  the  SE1  Hartree^-Fock  wave  function  {Sh  73)  derived  from 
the  Feshhach-Lomon  potential  increased  the  absolute  spectrum 
about  10%  for  b=1.70  fm  (b=1.70  fm  gave  the  best  agreement 
with  the  measured  40Ca  charge  form  factor  for  g<1  fm~l), 
however  the  relative  spectrum  was  reduced  less  than  3%.  This 
bears  out  the  insensiti vity  of  the  relative  spectrum  to  the 
choice  of  nuclear  wave  functions. 

Two  improvements  in  the  GD5  model  -  a  more  consistent 
choice  of  km  and  V7,  and  the  use  of  new  40Ca  photoa t sorption 
data  -  resulted  in  our  finding  that  the  GDR  model  relative 
photon  spectrum  is  comparable  in  magnitude  to  the  CHO  model 
spectrum  within  U%  for  55<k<85  Mev.  This  result  further 
indicates  that  the  relative  spectrum  is  to  a  large  extent 
independent  of  the  nuclear  wave  function.  The  GDB  model  thus 
predicts  absolute  and  relative  capture  rates  which  are  too 
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high  by  up  to  a  factor  of  two,  however  both  ordinary  and 
radiative  capture  rates  are  in  substantially  better 
agreement  with  experiment  than  CHO  model  absolute  rates.  The 
discrepancy  which  existed  prior  to  1977  between  experimental 
and  theoretical  determinations  of  the  4°Ca  photon  asymmetry 
has  largely  been  resolved  by  a  recent  SEEL  experiment 
(Ha  77)  .  Our  calculation  of  the  CHO  model  asymmetry  is  in 
agreement  with  previous  calculations  (Eo  65,  Eo  74)  and  with 
this  latest  measurement.  The  GDE  model  photon  asymmetry  for 
4°Ca,  a  quantity  not  previously  calculated,  was  found  to  be 
nearly  identical  to  the  CHO  model  asymmetry. 

The  modified  closure  calculation  presented  in  Chapter 
VI  indicates  a  possible  means  of  fixing  km  by  appealing  to 
experimental  phot cabsorption  cross  sections.  Though  our 
calculation  contains  a  number  of  approximations,  the 
underlying  ideas  can  be  pursued  further  if  one  is  interested 
in  improved  accuracy.  An  undoubtedly  temporary  limitation  in 
this  regard  is  imposed  by  contemporary  uncertainties  in 
measured  photoabsorption  cross  sections  used  to  determine 
the  modifying  factor  (1  +  S).  Eesults  obtained  for  the 
ordinary  capture  rate  by  modifying  the  classical  dipole  sum 
are  encouraging,  and  focus  attention  on  the  need  for 
considering  exchange  currents  in  both  ordinary  and  radiative 
muon  capture. 

As  previously  mentioned  in  our  discussion  of  motivation 
for  the  GDE  model,  nucleon- nucleon  correlations  play  a  major 
role  in  determining  the  position  of  the  giant  dipole 
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resonance,  and  are  capable  of  producing  substantial  changes 
in  EMC  observables.  The  shell  model  nuclear  matrix  elements 
described  in  Appendix  D  can  be  decomposed  into  one  and 
two-body  parts.  As  it  turns  out,  the  value  of  a  matrix 
element  depends  to  a  large  extent  on  the  partial 
cancellation  between  the  one  and  two— body  parts.  In  a  theory 
without  two— tody  correlations  such  as  ours,  the  degree  of 
cancellation  is  mcstly  determined  by  the  normalization  of 
the  wave  functions.  The  introduction  of  two-oody 
correlations  into  the  theory  might  possibly  affect  the 
cancellation,  thus  bringing  about  a  change  in  the  value  of  a 
nuclear  matrix  element.  An  investigation  of  two-body 
correlations  in  ordinary  muon  capture  by  McCarthy  and  Walker 
(Me  75)  has  shown  that  the  capture  rate  for  160  is  reduced 
by  30%  when  correlations  are  included- 

Lastly,  to  mention  a  closely  related  matter,  the 
possibility  of  nuclear  renormalization  of  the  hadronic  weak 
form  factors  has  only  recently  received  theoretical 
attention,  and  experimental  information  on  this  topic  is 
still  fairly  sparse  (Mu  76) .  The  approach  that  has  been 
followed  so  far  consists  of  attempting  to  retain  the 
framework  of  the  impulse  approximation  and  to  replace  the 
free  weak  couplings  with  effective  ones.  Indications  that 
the  axial  vector  strength  is  quenched  up  to  25%  for  medium 
and  heavy  mass  nuclei  (Br  78,  Os  79,  To  79)  have  strong 
implications  for  muen  capture.  Our  results  suggest  that  for  gA 
maximally  guenched,  both  ordinary  and  radiative  rates  are 
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substantially  reduced.  This  change  is  in  the  right  direction 
to  improve  agreement  with  experimentally  measured  rates. 
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Appendix  A 


Tree  Particle  Propagators  for  Maon  and  Hucleon 

for  the  muon,  the  momentum  space  representation  of  the 
relativistic  free  propagator  : 


So  C/t)  ^  77U 

M-z-  A-  1 

was  used  directly  in  the  evaluation  of  the  lepton  covariants 


In  the  case  of  the  nucleon,  a  free  propagator 
consistent  with  our  non-r eiat ivistic  Hamiltonian  was 
required.  Accordingly  it  should  satisfy  the  coordinate  space 
equation  : 


A-  2 


reduction  to  0(1/m2)  of  the  relativistic  Hamiltonian  for  a 
free  particle.  Since  the  initial  and  final  states  of  the 
nucleon  involved  in  the  EMC  process  are  energy  eigenstates, 
and  the  terms  comprising  the  effective  interaction 
Hamiltonian  have  position  and  time  dependence  : 


Hz  (x)  ~  Jtxp  (zA-x) 


Hw(x)  ~  JLXp  (KV'/^-x) 


A-3 


cur  problem  is  reduced  to  obtaining  a  suitable 
time-independent  propagator  for  the  nucleon  which 
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satisfies  : 

(£,r)„.sp-saf) 

A-4 

with  T—  E— in,  E  being  the  total  (relativistic)  energy  of  the 
nucleon  in  the  intermediate  state.  Thus  if  for  example  the 
neutron  emits  the  photon,  T=ly  =E^ +  k-m,  whereas  if  the  proton 
radiates,  T=Tp=E^-k- m.  Recognizing  that  Eg.  A-4  is  exactly 
the  equation  satisfied  by  the  usual  non-re  la  ti  vistic 
propagator  for  the  free  Schrodinger  equation,  one  can 
immediately  write  the  solution  : 


(*-p 


=  -L  /  dP^  exp  ( i  P-  (x-y)) 

Pz/2mT  + 


TJ  (*&  /-  Tyznn  A_5 

further  if  k  is  large  and  p2/2mT  is  small  compared  to  1,  the 
integrand  in  Eg.  A-5  can  be  expanded,  yielding  : 


<?.  (X-  v)  *  -L  {/  -If  +  V4  -■■■/  S(x-v) 
J  T  l  ZmT  (zmT)2-  J  J 


A-6 


Considering  the  neutron  radiating  diagram  specif ically ,  the 
second  order  S-matrix  takes  the  form  : 

’h  A  **  *&■ *  f'S.  'MS' ' 1 

x  S  (pz)j  Hw(z)  9i  (z)  A_  7 


ferforming  an  integration  by  parts  and  making  use  of  Egs. 


A-3  yields  : 
(2) 


$  */4 

+  (P-9- a-) 4 +...]}<&  (y) 

( ZmTN )x 

Similarly  one  obtains  for  the  proton  radiating  diagram  : 


A-8 
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-J dy  ty(y)  //w(y)  He  [/+ 

+  (P-k)?  +...JJ  9i(y) 

UmT,)z  J  7 

Note  that  the  correct  contact  limit  for  the  propagator 

(enclosed  in  braces  in  Eqs.  A  —  8  and  A~9)  is  obtained  by 

letting  the  nucleon  mass  go  to  infinity  : 


2ml} 


A-9 


Xv/rv  &o 

77l->  CO 


neutron  rout  cates 
pro  to  n  r  a.  diet  tes 


A-10 


At  this  point  the  initial  nucleon  was  put  on  mass 
shell ,  an  approximation  consistent  with  the  omission  of 
external  nuclear  fields  : 


-*  2 

m  ~  EL 
2m 


A-11 


Expanding  *lp  and  in  the  variable  p2/2mk  and  keeping  terms 
formally  to  0(1/m)»  the  free  propagator  is  approximated  : 


_/ _ 

(ny-v) 


+  -  x  AiiL  z 

2m  (77pr?)z  m  m  ( npr?) 


k'P 

k  m  mk 


proton 

radiates 


neutron 

radiates 


A-12 


As  it  turns  out,  higher  order  terms  in  the  expansion  of  the 
propagator  beyond  the  first  P/m  term  are  not  required  if  the 
squared  matrix  element  is  kept  to  0(1/m2)  ;  in  fact  keeping 
terms  in  the  squared  matrix  element  to  second  order  meant 
that  in  most  cases  only  the  leading  contact  term  in  the 
expansion  of  the  propagator  was  actually  used- 


Appendix  B 

BMC  Effective  Hamiltonian 


Cur  effective  Hamiltonian  has  the  overall  structure  . 

H tU  =  /y  O-  0L-  t)  771 fa 

with  : 


B-1 


7?l  *  77[,  +  01' 77[z 


3-2 


lo  facilitate  manipulation  by  the  algebraic  squaring  routines, 
77£  was  arranged  into  the  form  : 

£  ?-(/){£ tZ  (£i  tPj.  +&)&■&}. 

L  /*/  J  Li*/  i*NH  J  B-3 


where  F.  and  are  functions  of  the  vectors  6,  k,  y>,  P/m,  and 
O' .  The  coefficients  E.  ,  D*  ,  and  C :  correspond  to  terms 
of  order  1,  1/m,  and  1/m 2  respectively.  Table  III  lists  the 
terms  appearing  in  77£,  r  while  table  IV  lists  the  terms  in  ^2- 
A  •' Y **  appended  to  a  coefficient  name  denoted  that  the  factor  y 

A  4 

M-V  appeared  in  the  coefficient.  These  "Y”  coefficients 
were  segregated  because  an  integration  on  y  had  to  be  done 
after  the  algebraic  squaring.  A  complete  list  of  the 
coefficients,  as  used  by  the  numerical  program  NEKBMC, 
appears  below.  The  FORTRAN  symbols  are  explained  in  table  V. 
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i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

13 

14 

15 

16 

17 


Table  III 
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Terms  in 


Non-Zero  Coefficients 


*  A 


u,  e-  v  *  cr 

D1 

Cl 

A  ->  -* 

A  t-P  x  a/m 

C2 

X  6  -  P  X  W/77 

C3 

A 

e- cr 

E4  D4 

C4,C4Y 

e-p/m 

D5 

C5 

A  A 

e-  v 

C6 

A-P  i-7U 

C7 

P-  P  £■  t/fif 

C8 

P-  V  £•  7/fn 

C9 

A  A  A 

v-cr  e-y 

D1 0,D1 0Y 

Cl  0 ,  Cl  0Y 

7  A  A 

&■  cr  6-V 

Dll  ,D11  Y 

Cll ,C1 1 Y 

P-  7  £■  v/m 

Cl  3 ,  Cl  3Y 

P-7  P-e  /mz 

Cl  4 

P-e  v-  7/m 

Cl  5 

P-e  A- 7/m 

Cl  6 

P-e  P-M/m2- 

Cl  7 

Combinations  of 

vectors  appearing  ir 

i  th  e 

expression 

for  777 i ,  with  the  effective^  Hamiltonian  given  by 

=  &(/-£.’ v)(7rij  +  (7i' m*-) fa 

The  non-vanishing  terms  cf  order  1 ,  1/m,  and  1/m2 

are  denoted  by  the  E,  D,  and  C  coefficients 
respectively- 


i 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

37 

38 

39 


Table  IV 


Terms  in 
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G  • 

l 

Non-Zero  Coefficients 

*  r. 

X  €  X  CT 

E20  D20 

C20,C20Y 

•  A  1 
z  e  x  y 

C21 

z  A- P  ex  o' /m 

C22 

X  P-7  ex 

C23 

^  v7/7  e  xtr/sn 

C24 

xd  v-axp/m 

C25 

ze  A-  o' x  p/m 

C26 

x  R  V-  or  xe 

C27 

Z  A  €•  X  P/77! 

C28 

•  A  A  77  / 

ZC’V  P X  (7V777 

C29 

.  A  A  a  _> 

z  £•  v  V  x  cr 

C30 

z  e-  v  A  x  cr 

C31 

z  P-  cr  i  x  v/m 

▼ 

C32 

z  P-  &  Px  e  /m2 

C33 

z  v-vPxe/771 

C34 

z  A- &  fixe/m 

C35 

z  e  x  P/777 

D37 

C37 

z  P’€  v  x  cr / 777 

C38 

A 

e 

E39  D39 

C39,C39Y 

Combinations  of 

vectors  appearing 

in  the 

expression  for 

m4 with  the  effective  Had. It onian  given  by 
Htf*=  /Ui-  <£•  v)(7n<  +  ai;  z 

The  non- vanishing  terms  cf  order  1 ,  1/m,  and  1/m2 

are  denoted  by  the  E,  D,  and  C  coefficients 
respectively. 


1 

40 

41 

42 

43 

44 

45 

46 

47 

48 

50 

51 

52 

53 


Table  IV  (continued) 
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Terms  in 


G.j  Non-Zero  Coefficients 


£  P-  d/m 

D40 

C40  ,C40Y 

k  e-a 

D41 

C41 

-»■  A  A 

O’  6-Y 

D42 

C42 

A  ->  A 

e  c- v 

D43,D43Y 

C43,C43Y 

p  S  •  r/m 

C44 

e  p-  k  / 771 

C45 

e  P-  P/m2- 

C46 

i  p.  ?/m 

C47 

a  e-p/m 

D48 

C48 

A  £  •  P/m 

C50 

PG-  V /m 

C51 

P  £  ■  P/mz 

C52 

£  O'- A 

D53.D53Y 

C53,C53Y 

Table  V 


lable 


FORTRAN  Symbol  Dictionary 


FORTRAN  Symbol 


Variable 


x 
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LPL 

X* 

LM 

X- 

MMU 

vy 

M 

7T1 

MPI 

771k 

MUP 

Xp 

MUN 

Xh 

MUA 

Xp  Xn 

MUS 

Xp  'K* 

MU 

I+Xp-Xh 

KS 

/A/ 

K 

A 

k 

NUS 

/V/ 

NU 

A 

V 

FSDK 

A  A 

$  term  flag 

FVEL 

velocity  term  flag 

FMR 

muon  radiating  diagram  term  flag 

KM 

Mm 

FCP 

contact  propagator  term  flag 

GM 

S* 

GP 

Sp 

Dictionary  of  FGETEAN  symbols  appearing  in  the 
computer  listings  of  the  sguared  matrix  element- 

Table  V  (continued) 


FORTRAN  Symbol 

Dictionary 

FORTRAN  Symbol 

Variable 
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km /2m 
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(2m) 

RNO 

(rn^'V) 

RN1 

(Vz/Zm)  / (77u-v)Z 

FNM 

=0  (unused  flag) 

appended  G 

denotes  complex  conjugation 

F4  =  -GA*|_Pl*FMN 
F20=GA  Al.PL  *F  MM 

E39=GS* ( I . -L  Pi  *FMR ) -GV  *L PL  *FMR 
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Appendix  C 

B2!C  Matrix  Element  Squared 


Ihe  REDUCE2  programs  SQUARE1,  SQUARE2,  and  SQUARE3  performed 
the  algebraic  squaring  of  ^  •  with  output  from  SQUARE3  being 
done  in  the  FORTRAN  language.  The  complete  expression 
for  was  not  squared  all  at  once,  rather  it  was  divided 
into  segments  according  to  powers  of  1/m  (via  the  E,  D,  and 
C  coefficients  of  Appendix  B) ,  and  the  indivdual  segments 
processed.  Listed  below  are  the  parts  of  2  which  were 
calculated  : 

1.  M  SO  1  :  (0  ( 1 )  )  2  +•  0(1)  xO  (1/m) 

2-  MS012:  (0  ( 1) +0  (1/m)  )  2 

3.  MS2  :  0  (1)  xO  (1/m2) 

4.  MS 3  :  0  ( 1/m)  xO  ( 1/m2) 

The  output  files  MS01,  MS012,  M S2 ,  and  MS 3  are  reproduced 
below.  Products  of  nuclear  matrix  elements  are  denoted  : 


C-1 


where  the  expressions  POLn  are  defined  in  Appendix  D. 
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Appendix  D 


Nuclear  Matrix  Elements  in  the  Closure  -  Harmonic  Oscillator 

Model 

In  the  standard  theory  of  radiative  muon  capture  the 
evaluation  of  : 

<G/><02>  =  Z  (bill  T-  (j)  txp  (-IS'Aj)  0t{j) 

s 1 

X  <b/Z  T-  (k)  exp  (-iS  Ak)  OzM/^f  ^ 

is  a  principal  concern,  where  |a>,  |b>  are  the  initial  and 
final  nuclear  states,  respectively,  and  0 ( j)  is  an  operator 
referring  to  the  j*th  nuclecn.  To  evaluate  products  of 
nuclear  matrix  elements  of  the  form  D-1  the  closure 
relation  : 

£  Ib'Xbl  =  1 

b  D-2 

was  utilized  and  further  it  was  assumed  the  initial  state  is 
a  spherical,  spin  zero  nucleus  which  can  be  described  by  the 
shell  model  as  a  set  of  completely  filled  shells. 

Accordingly,  a  Slater  determinant  was  chosen  for  i a>  : 

/«>  -  -U.Z (l)  9kz {&)...  9k, (Z) 

VA!  p  '  £>-3 

Here  A  is  the  nucleon  number  and  P  is  the  coordinate 
permutation  operator.  is  the  wave  function  for  a  particle 
in  the  shell  model  state  Hence  : 

KoXOz'y  =  zjZKZjotoJti-Z  <h/exp{-is-A)0,/X,y 

X  OizUxp  (~iS -A)  0z  / X/}  }  d-  4 
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uhere  the  index  [%z)  takes  on  all  values  of  the  coordinate 
space  quantum  numbers  describing  occupied  proton  (neutron) 
states. 

In  our  study  of  EMC  the  following  products  of 
expectation  values  arose  : 

1.  <1><1>* 

2.  A-<P><1># 

3.  <?>•<£>* 

4.  A- <P>  B*  <P> 

5.  <?-PXl/ 

6.  A»  <P  B •  P><  1> 

7.  AXPXB*"?  C*"p>^ 

8.  <P>-  <P  !•?>* 

S.  <  A  •  P> <~P  •  ~P>* 

where  P=-i7  and  A,  B,  C  are  vectors  not  depending  on  the 
nucleon  coordinates.  Expressions  for  the  combinations  1  and 
2  can  be  found  in  the  literature  (Lu  63),  whereas  the 
quantities  3  through  S  have  not  appeared  in  previous  EMC 
studies.  Expressions  for  the  latter,  which  appear  here 
because  we  keep  all  0(1/m2)  terms,  were  derived  as  follows. 
Harmonic  oscillator  functions  were  taken  for  the  single 
particle  wave  functions  and  the  evaluation  of  expectation 
values  was  done  in  Cartesian  coordinates.  The  unnormalized 

wave  functions  are  : 

^ny7iy7iz  ~  //n* (X) Hny (y)  f/nz (Z)  ^ 

with  Hn  a  Hermite  polynomial  of  degree  n.  Each  expectation 
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value  required  that  an  integral  of  the  form  : 


le  evaluated.  The  necessary  integrals  were  derived  from  the 


relation  : 

oo 


-oo 


D-7 


hy  successively  differentiating  under  the  integral  sign  with 
respect  to  s  to  generate  each  term  in  the  polynomial 

The  final  results  involved  polynomials  in  sb 
irultiplied  by  the  Gaussian  factor  e xp (- (sb)  2/4)  . 

The  products  of  expectation  values  containing  an  odd 
number  of  P  operators  {2,1, Q, 9)  were  evaluated  by  expressing 
them  as  combinations  of  other  expectation  values  given  in 
the  list  above.  A  parity  argument  relates  the  expectation 

values.  By  way  of  illustration  of  this  relation  the 

-►  *■ 

evaluation  of  Ee<?><1>  is  considered  below.  The  extension  of 
the  argument  to  the  evaluation  cf  the  products  listed  as  7, 

8,  and  9  is  straightforward. 


Re  <P><!>*~  i  {<?>*«>  +  <P><!>*} 
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D-8 

Note  that  there  is  do  one-body  contribution  to  Ee<T>< 1 >* 
since  the  parent  nucleus  is  assumed  to  be  spherically 
symmetric.  Now  change  r  to  — r  in  the  second  term,  assuming 
the  states  2/?2z  have  definite  parity  : 

/Re  <?></>* -  -Z  $  /<2,/&xp  (-itA)/A*>/2' 

(2z  / txp  Us- a)  /A/)  Ot/ txp  f-i  s-Z)  P  /2Z> 

-  <2/  / zxp  U  s-a)/Az>  Oz  Uxp  (-/  T  a)  P  /2i> 

D-9 

If  N=Z  and  states  2/  are  identical  with  states  2z ,  2/  and  2z 
can  be  interchanged,  making  the  last  two  terms  vanish  : 

/Re  <P></>*=  ~sZ  I  Oil  exp  (-1  s-/0  /2z  >/  * 

2/^2z 

~  TB  D-  10 

Ihus  Ee<?><1>  is  expressed  in  terms  of  the  two-body  part  of 

i<  1>  I  2« 

Ihe  computer  language  IEDUCE2  (He  73)  was  used  to  sum 
the  polynomials  involved  in  each  product  of  expectation 
values,  final  results  for  the  products  of  nuclear  matrix 
elements  are  : 

/).  /< v>/* =zo{/- a  +£-£+$ *&)} 

s  POU 

Z).  /Re  A-<P></>*  *  ~  /0A'2_(/+fy- 

~m  77i  o  ou  yfO 

2  -  A’JL  x 

2m 

{2-("pifw  * nL3 
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“  A'B*x  P0L4  -  A-  S  3* S  x  POLS 

2771  2771 


S).  <P-P>  </>*=  M 

77lz  M)z  ZO  Z40 

s  POL  6 


6).  <AP3-p) </>*  =  ZO  A- 


771 


■  z  A- 1 B-JL  ( /-2L  -f-Z.  2- B  A  P8 )exp(-?i/z) 
m  m  z  /6  76  ZS6  ' 


-  A-B  x  P0L7  ~  A-S  B-_S_  x  POL 8 

2772  2771 


7).  Re  <'AP>*S-PO-P>  =  Ap  {b-CxPOL7-ZOb2-  ~Pl.  C-_S_xP0L8} 

771 3  2771  (mb)Z  2771  2771 

y-  B-S  {I-CxPOlf  -20A'C  -A'l  carols} 

2m  (tti6)z  2771  2771 

+■  C-_L  lA-B  * P0L4  -Z0_±s  -A^s.  .SvS  *  POLSj 
Zm  (7716)*  Zm  Zm 

7  A-S  3-S  C-S  x poL Z 

8m3 


6).  Re  <p>-  <PA-py  =  A-l  l P0L7- 100- xz  ~Ju_ x POLS  -f- POL 3 
m3  2771  (mb)  4 771 

+  POLxf-  -B.XPOLS  P°AZ] 

4mz  2771 


j 


. 
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?)■  /Re  <AP><PP>*  =  A-S  { POL/  +  2xP0L4  -  /OO 


(mb)z 


( POLZ-  POLS)} 

2m z 


fO).  POL  (S  =  20 


D-11 


In  our  GDE  model  calculations,  the  elastic  form  factor 
**  and  the  unretarded  dipole  part  D  of  the  main  nuclear 
matrix  element  squared  were  also  required.  Calculation  of 
these  quantities  proceeded  as  above  for  the  matrix  elements, 
lhe  results  are  : 


J7T  i=/  /-/ 


D-  12 


D-1 3 


j 


Appendix  E 


Bain  nuclear  Matrix  Eleaent  Squared  Csing  the  SBL  Save 


Function 


Taking  the  ground  state  wave  function  to  be  a  Slater 
determinant  : 


fay  “  ZL  9ic,  CA,)  9bcz{Az)..,  (Aa) 

E-1 

and  substituting  the  SBL  (Sh  73)  single-particle  wave 
functions  : 

Pic  (a)*  ^KjCJiTL-t  (a)  *  ZL  ZL  Cm//  (Zmz  z  m*  limy  Rmj(a) 

J  '  M  TTtjTTls  J M  J 

X  Yt  *{J2a)  OCs  OCt 


=  Z  ZL  CnmjL;  <7 m*  z  Trts/jvy}  /M/rrii)  lzmZ?lz7nZ> 

M  77U ,7/U  7  v  7  E-2 

in  Eg.  D-1  yields,  in  closure  approximation  ; 


</><»*-  z  a///to-raA)<J.rty 

A/  A/  ?Az 


E-3 


where  Af  (iz)  takes  on  all  values  of  the  quantum  numbers  (Nljmj) 
of  states  occupied  by  protons  (neutrons).  The 
tNMij  are  constant  coefficients  determined  by  Shao  et  al 
(Sh  73)  and  reproduced  in  table  VI ,  (^771/  !j  77ljy  are 

Clebsch-Gordon  coef  ficients,  and  ,  H  771$) t  Hl  77tfZ  are 

space,  spin,  isospin  wave  functions,  respectively.  Working 
in  spherical  coordinates,  the  angular  momentum  algebra  can 
be  handled  using  standard  techniques  (Ed  57),  so  that  a 
straightforward  calculation  yields  : 


186 


J 


Table  VI 


Oscillator  Function  Coefficients  .  for  ^Ca 

(Normal ization=l ) 


Nlj  IS  1/2 

M 

2S  1/2 

IP  3/2 

IP  1/2 

ID  5/2 

ID  3/2 

1 

0.9575 

0.1709 

0.9867 

0.9994 

0.9947 

0.9951 

2 

-0.2634 

0.8223 

-0.1247 

-0.0030 

-0.0730 

-0.0692 

3 

-0.1141 

-0.5005 

-0.0958 

-0.0289 

-0.0527 

0.0684 

4 

-0.0239 

0.1592 

-0.0419 

-0.0163 

-0.0500 

-0.0160 

5 

0.0059 

-0.1262 

0.0000 

-0.0066 

-0.0083 

0.0119 

6 

0.0086 

0.0536 

Table  VI-  Oscillator  function  coefficients  tot  the 

40Ca  wave  function  of  Shao  et  al  (Sh  73)  .  The 
single  particle  wave  functions  are  normalized  to 

one. 
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</></>  sZ'ZMz  Z.  Z  Cu.M.ty,  C/bMMJUft  Chm'Ij,  CtoMz'Ay,  (ZZ+l)  (Zf,  H) 
W,M*  N*^j*  X  (2j*+i)  ItolMihi  Z*u,mL/u  (Jj_  ^q)Z  e-4 

The  last  factor  is  the  square  of  a  Kigner  3-j  symbol  and  the 
radial  integrals  are  represented  as  : 


I  Ml  M  'I'L  *  J '  R  Ml  fa)  R+M'l’fa) JL  (SA.)  A  Z<A 


E-5 


with  (s r)  a  spherical  Bessel  function-  Bor  the  harmonic 

oscillator  potential.  Eg.  4.97  of  DeForest  et  al  (De  66) 


cives 


Imim'i'l 


(uh)H  [z 


Z  L 

/?i\z 


z  J.  M 

exp  (-  %)(( 'M-/)! Ym'-/)/)  1(r(M+t+i)  r(M  ti  'H)Y 


M-!  N+N' 

xz  Z 


/ 


n=o  n'zo  H!N'!  (M-N- OKm-N-I)!  r(N+ Z+3/Z)  T(N'-t Z'f  %■) 

x  fCzCl-Z-I'-zn-zn'))  l+%)  *1%.) 

E-6 

where  7£  -  (sb)  ,  T  is  the  gamma  function,  and  F  is  the 
confluent  hyper ceometric  function. 

Taking  <ZNfAlj  -S//m  recovered  the  simple  wave  functions 
used  in  Appendix  D  and  allowed  us  to  check  that  for  this 
simpler  case.  Eg.  E-4  reduced  to  the  appropriate  expression 
which  appears  as  Eg.  4.10  in  Luyten  et  al  {Lu  63)  • 

Inspection  of  Fg.  E— 6  shows  that  the  evaluation  of  rg. 
E-4  must  lead  to  a  result  of  the  form  : 


2 

z-z  Cz.  +a,  if  +  yV-3  exp  (-  V/z) 


E-7 


hence  our  task  was  to  calculate  the  a^-  -  This  was 
accomplished  by  a  FCETEAN  program  which  was  checked  insofar 
as  possible  ty  letting  •  thu£  recovering  the 
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corresponding  expression  obtained  for  the  simple  wave 
functions.  Our  final  result  was  : 


</></>  =  10- (zo- 3.0571  +  SJ5 if- U3 71**0. 4/4 q  )  e 


E-8 


For  the  GDR  mcdel#  it  was  necessary  to  also  calculate 
the  unretarded  dipole  part  D  of  the  main  nuclear  matrix 
element  and  the  elastic  form  factor  F^.  These  calculations 
proceeded  in  a  wholly  analogous  fashion  to  the  one  above  and 
hence  only  the  final  results  appear  here. 


D  =  /3.0S  7?Z  E_9 

2 

Fa  =  (/-  0. 346  if+  0.  0366  if-  0.  004/6  e-  1  o 


Appendix  I 

A  Consent  oje  Fearing *s  Thsocei 

The  theorem  proved  by  Fearing  {Fe  75)  states  that  : 

•'for  the  standard  theory  of  radiative  muon  capture 
including  all  weak  couplings  except  gs  and  all  usual 
diagrams,  the  photon  asymmetry  OC  and  the  circular  polarization 
ft  satisfy  <**=+ 1+0  ( 1/m2)  and  p  =+  1+o  ( 1/m2)  , 
where  m  is  the  nucleon  mass." 

Ihe  standard  theory  referred  to  (Eo  65,  Fe  66)  considers  the 
muon  and  nucleon  to  propagate  freely  in  their  intermediate 
states.  In  what  follows  the  question  of  whether  the  theorem 
remains  valid  when  external  potentials  are  incorporated  in 
the  theory  is  examined. 

Consider  the  Feynman  diagram  in  figure  1(a).  The 
corresponding  amplitude  can  be  written  ; 

Ai/  ~  y  K  ( ¥5)  So  T 

F-1 

where  S0  is  the  free  muon  propagator  and  T*  the  hadronic 
current.  This  amplitude  forms  part  of  the  standard  theory 
considered  by  Fearing  (Fe  75).  If  figure  1  (a)  is  altered  so 
that  the  intermediate  state  muon  interacts  with  some 
external  field,  then  to  first  order  in  the  coupling  to  the 
external  field  the  new  amplitude  is  : 

M,  -  ft  £  a-fs)£S'  (f+/*-k)Hr(fUS.(/i-k)  Ht)  ^ 

with  T(q)  the  external  interaction  vertex.  Note  that  to  all 
orders  of  a  lippmann-Schwinger  expansion  of  the  muon 
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propagator,  the  rightmost  factor  will  he  Se  (This  fact 

will  be  used  later  when  our  argument  is  generalized.)  All 
other  diagrams  of  interest  (figs.  1(b) -1(e),  modified  to 
include  the  external  field)  may  be  lumped  into  the 
amplitude  : 


Mz  ~  Yk  (/-  /s)  Tfa 


E-3 


/A 


^^^■6  ^  represents  photon  emission  by  the  hadronic 
current.  Using  Egs.  E-2  and  F-3  yields  : 

M2-Z  /M,+Mz/z 

spins 

-  Tr  Hf> (his)  (fa  (at/L-k) rpSc Cu-A)Jt  TWy?  UtnuYMXr’fl 

x  (ffapJdms.  to+/i-k\jT+R%')  Is  (/-Ys)jL 

^  F  —  4 


laking  the  muon  at  rest  and  using  the  explicit  form  for  the 
free  muon  propagator  : 


60  (u-fe)  = 

r  (/c-k)z  -  7?JuZ 


-2k  77^ 


F-5 


yields,  after  some  algebra 


(his)  (fa  Tcpl  M/M)  ) 

y  fl+  ft. ■ .?) tin. ■  £ *(X + Is) [r*p S' (?  +/V-M)  T  +R  l\)Y(  (hjs) 


2t7J~ 


F-6 


Now  consider  the  factor  (S0  (g^_k) /’(q)  )  •  If  higher  order 
terms  from  the  expansion  of  S  were  kept  in  M/  ,  this  factor 
would  be  different.  Independent  of  its  explicit  form 
however,  it  can  always  be  decoirposed  into  a  sum  E+C  such 


that  : 
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Yr  £  =  £  Ys 
YsO  *  -  OYs 

F- 7 

representing  the  factor  in  question  by  the  sum  E+0  and  using 
the  commutation  rules  of  Eg.  F  —  7  allows  us  to  write  ; 

M  - r^rf(/-)fs)(££lt-02^J(j±M  &■£  Tf+R%)(/tK-s) 

77J*  77^  Z  r 

x (oi-  e'f2t  £*- Jr  omJ  r‘sf-Rs% )  it 

77ju.  77)*,  6  2777?  F-8 

Now  look  at  the  first  term,  which  is  the  square  of  the  muon 
radiating  amplitude.  Taking  explicit  representations  for  the 

A  4  4 

gamma  matrices  and  using  6  =  (i-iXj) //2 ,  the  term  becomes  : 

ALr  =  Tr£  Y*  (/-Ys)  £(J±M  T^2+  (/-  £Lk)Of  S-k)  £*TS JtL  A 

7  z  2  m,  7. tf» 

tTrYYp  (/-  Ys)  0  (/+Y.)  T?X-  Of  oL-k)(h  Ik)  cffs  Ys 

2  2777v?7)Z  F"S 

Contact  is  made  with  the  work  of  Fearing  (Fe  75)  ,  in  which 

* 

the  muon  was  taken  to  propagate  freely,  by  setting  E=E  =1 
and  0=0  =0.  In  this  case  M;  has  the  overall  f  actor  A^{  1* 5-k) 
and  thus  has  X- 1.  In  the  general  case  where  the  muon 
propagates  in  some  external  field  however,  0  will  not 
vanish.  For  this  general  case  then,  the  muon  radiating 
diagram  will  contribute  to  Xi-  1.  Thus  the  conclusion  is 
reached  that  Fearing* s  theorem  does  not  apply  to  theories  of 
radiative  muon  capture  which  include  propagation  of  the  muon 
in  external  potentials. 


